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Abstract 

Background: The Fox gene family is a large family of transcription factors that arose early in organismal evolution 
dating back to at least the common ancestor of metazoans and fungi. They are key components of many gene 
regulatory networks essential for embryonic development. Although much is known about the role of Fox genes 
during vertebrate development, comprehensive comparative studies outside vertebrates are sparse. We have 
characterized the Fox transcription factor gene family from the genome of the enteropneust hemichordate 
Saccoglossus kowalevskii, including phylogenetic analysis, genomic organization, and expression analysis during early 
development. Hemichordates are a sister group to echinoderms, closely related to chordates and are a key group 
for tracing the evolution of gene regulatory mechanisms likely to have been important in the diversification of the 
deuterostome phyla. 

Results: Of the 22 Fox gene families that were likely present in the last common ancestor of all deuterostomes, 
5. kowalevskii has a single ortholog of each group except FoxH, which we were unable to detect, and FoxQ2, which 
has three paralogs. A phylogenetic analysis of the FoxQ2 family identified an ancestral duplication in the FoxQ2 
lineage at the base of the bilaterians. The expression analyses of all 23 Fox genes of 5. kowalevskii provide insights 
into the evolution of components of the regulatory networks for the development of pharyngeal gill slits (foxC, 
foxLl, and fox/), mesoderm patterning {foxD, foxF, foxG), hindgut development (foxD, fox/), cilia formation (foxJi), and 
patterning of the embryonic apical territory (foxQ2). 

Conclusions: Comparisons of our results with data from echinoderms, chordates, and other bilaterians help to 
develop hypotheses about the developmental roles of Fox genes that likely characterized ancestral deuterostomes 
and bilaterians, and more recent clade-specific innovations. 

Keywords: Hemichordate, Saccoglossus kowalevskii, Fox cluster, Deuterostome evolution, Gene regulatory networks, 
Gill slits, Gut patterning, EH-l-like motif, FoxQ2, Fox genes, Forkhead, Fork head 



Background 

The development of animal body plans and associated mor- 
phological innovations are a result of genetic and cellular 
mechanisms acting in space and time. Developmental regu- 
lation of these mechanisms has many layers of complexity, 
and involves interacting suites of transcription factors that 
form core conserved regulatory kernels [1]. Analyzing these 
transcription factors and their genetic interactions is there- 
fore essential to understand development, and comparative 
studies between species can help us to understand the 
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evolutionary traits of developmental programs/networks 
and how they arose during evolution [2-6]. 

One large transcription factor family with key regula- 
tory roles is the Fox (Forkhead box) transcription factor 
family. Fox genes encode transcription factors contain- 
ing a fork head helix-turn-helix DNA binding domain of 
100 amino acids [7-12]. The conserved protein sequence 
encoding the DNA binding domain was described in 
1990 by comparative analysis of the Drosophila melano- 
gaster 'Forkhead' protein [13] with the HNF-3 protein 
isolated from rats [14,15] by Weigel et al. [16]. In the 
20 years since their discovery, a large number of Fox 
genes have been characterized in a phylogenetically broad 
range of animals, including choanoflagellates, yeast, and 
fungi (reviewed in Larroux et al. [17]) and a unified 
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nomenclature of 15 Fox families (alternatively named clas- 
ses or subclasses) was established in 2000 [18]. Five more 
families have since been added: FoxAB [19-21], FoxQ 
[22-24], FoxP [25,26], and the vertebrate specific groups 
FoxR [27,28] and FoxS [29]. Four of these families have 
subsequently been further divided: FoxL into FoxLl and 
FoxL2, FoxN into FoxNl/4 and FoxN2/3, FoxQ into FoxQl 
and FoxQ2, and FoxJ into Foxjl and FoxJ2 [30]. This has 
led to the identification of a total number of 24 Fox families, 
making it possible to compare their orthologs in different 
species to gain insights into the evolution of this large tran- 
scription factor family and their roles in metazoan develop- 
mental programs. 

Fox genes probably arose by serial duplication from a 
single ancestral Fox gene present in the fungal/metazoan 
ancestor (stem opisthokont) [17,18]. This evolutionary 
history is reflected in the clustered arrangement of some 
of the Fox genes in animal genomes [31-33]. By com- 
parative genomic analysis two Fox gene clusters have 
been proposed to be present in stem bilaterians; a foxD- 
foxE cluster and a foxLl-foxC-foxF-foxQl cluster [31,33]. 
The latter is of special interest since its conserved link- 
age may be correlated with its function in mesoderm de- 
velopment [31,32,34]. Much of the literature on Fox 
genes focuses on medically relevant developmental roles 
using data from a narrow range of vertebrate model sys- 
tems including only a few invertebrates like Drosophila 
melanogaster [35-38] and Caenorhabditis elegans [39-42]. 
Recently, a more extensive evolutionary comparative lit- 
erature has begun to emerge; new data from animals such 
as elasmobranchs (dogfish) [34], echinoderms [20,43-48], 
cephalochordates [21,49-55], urochordates [56-60], lopho- 
trochozoans [32,61], cnidarians [19,62-65], and sponges 
[66] make it now possible to investigate the deeper evolu- 
tionary history of this transcription factor family. 

In this study, we have characterized the full Fox gene 
complement of the enteropneust hemichordate Saccoglossus 
kowalevskii to contribute to this discussion. Hemichordates 
are a deuterostome phylum, sister group to echinoderms, 
which together form the Ambulacraria [67-70]. Hemichor- 
dates share many organizational features with chordates 
such as a bilateral body plan with a conserved anterior pos- 
terior patterning gene regulatory network [71-73]. Their 
anterior gut is perforated by pharyngeal gill slits, likely 
homologs to those of chordates [74-79] and they have a 
nervous system with both diffuse and central organizational 
elements [71,80-87]. 

In our study, we identified 23 Fox genes in S. 
kowalevskii and analyzed their phylogenetic relation- 
ships, genomic cluster organization, and spatiotemporal 
expression patterns during early development. The ex- 
pression analysis of all 23 genes gives insights into the 
evolution of components of the regulatory networks for 
pharyngeal gill slits (foxC, foxLl, and foxl), mesoderm 



patterning (foxD,foxF, foxG), hindgut development (foxD, 
foxl), cilia formation (foxjl), and patterning of the apical 
territory (foxQ2). 

Methods 

Embryo collection 

Adult S. kowalevskii were collected at Waquoit Bay, 
Massachusetts in September. Oocyte ovulation and 
fertilization were carried out as described previously 
[88]. Embryos were staged by the normal tables of 
Bateson [74,89,90] and Colwin and Colwin [91]. Em- 
bryos were cultured at 20°C. 

Identification and cloning of Fox genes 

Numerous Fox genes were identified by screening expressed 
sequence tags (EST) [92], and expression patterns of 
foxQ2-l, foxG (bfl), foxA, and foxC in select develop- 
mental stages have previously been published [71,79,93]. 
To identify and clone additional genes, we screened the 
5. kowalevskii genome-trace archive at NCBI and our 
arrayed EST clone libraries [71] by bidirectional blast. 
Genes not in our EST libraries were cloned by PCR 
from cDNA prepared from a variety of developmental 
stages, using RNAeasy (Qiagen) for RNA extraction and 
Superscript III (Invitrogen) for cDNA synthesis, and cloned 
into the pGemT easy vector system (Promega). Primers 
used to clone partial fragment of FoxJ2/3: 5'-CAATG 
GACTGGCTGCCACAACTA-3 ', 5 ' -GTGTGAAGAACT 
GATTGAGTGAATTTGC-3 ' . 

In situ hybridization 

In situ hybridization was carried out as described in Lowe 
et al. [88] with the following modifications: Proteinase K 
treatment was carried out at 10 ug/mL for 5 min at 37°C. 
Acetic anhydride treatment at 250 uM for 5 min at room 
temperature (RT) followed by a 500 uM treatment for 
5 min at RT. 

Sequence retrieval 

Reference sequences for the alignment were assembled 
from a variety of metazoans; cephalochordates [21,94], 
sea urchin [20], cnidarians [17,19,62,63,95,96], and sponge 
[17]. For a list of GenBank accession numbers, see 
Additional file 1: Table SI, Additional file 2: Table S2 
and Additional file 3: Table SI. Sequences were aligned 
using ClustalW (EMBO-EBI). Trees were constructed 
in FigTree vl.2.3 [97] and further modified in Adobe 
Illustrator CS3. 

Additional software 

The DNASTAR Lasergene software package was used 
for sequence management, genome walking, and initial 
alignment. 
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Molecular phylogenetic analyses 

All genes in this study were assigned orthology by phylo- 
genetic analysis. Two types of analysis were carried out: 

(1) Bayesian analysis using MrBayes (v3.0B4) [98,99]; and 

(2) maximum likelihood analysis using the web-based 
PhyML server at Lirmm [100,101]. Alignment of the 
Fox (fork head box) domains was performed using 
ClustalW2 via the EMBL-EBI homepage [102]. 

Phylogenetic analysis of S. kowalevskii Fox proteins 

Bayesian analysis (MrBayes (v3.0B4), [98,99]) was carried 
out using the mixed amino acid substitution model apply- 
ing four independent simultaneous Metropolis-coupled 
Markov Chains Monte Carlo in two independent simultan- 
eous runs. N chains was set to 16, and the tree was calcu- 
lated on a 32 CPU cluster. The likelihood model was set to 
gamma rates = 4. A tree was sampled every 6,000 genera- 
tions for 53 million generations. The first 25% of the sam- 
pled trees were excluded via 'burnin' prior to consensus 
tree calculation. Saccharomyces cerevisiae Foxl was used as 
an outgroup. The trees converged to a standard deviation 
of 0.0109. Because of the size of the dataset, the maximum 
likelihood analysis was performed using the Approximate 
Likelihood-Ratio Test (aLRT) using the SH-like model 
[101]. The input alignment is comprised of 200 sequences 
with 88 characters (see Additional file 4: Table S4). 

Phylogenetic analysis of the FoxQ2 family 

The fork head box of FoxQ2 proteins from various phyla 
was used for this analysis, including sequences from the 
genome of the mollusc Lottia gigantea for which we 
identified four putative FoxQ2 genes by bidirectional 
blast and named foxQ2-l to foxQ2-4 (see Additional file 
5: Table S5). We only included proteins in our analysis 
which were short branching, and had an identifiable EH- 
I-like motif on the N-terminal or C-terminal side of the 
protein (see below) (Additional file 5: Table S5). The 
EH-I-like motif is outside the fork head box and was 
thus not part of the alignment. Bayesian analysis was 
performed using the Jones amino acid substitution model, 
applying four independent simultaneous Metropolis-coupled 
Markov Chains Monte Carlo in two independent simul- 
taneous runs. The likelihood model was set to gamma 
categories = 4 and gamma rates = invgamma. A tree was 
sampled every 500 generations for 1 million generations. 
The first 25% of the sampled trees were excluded via 
'burnin' prior to consensus tree calculation. Nematos- 
tella vectensis foxO was used as an outgroup. The trees 
converged to a standard deviation of 0.023. Maximum 
likelihood analysis was performed using the aLRT (SH- 
like model) [101]. The input alignment is comprised of 
28 sequences with 136 characters (see Additional file 6: 
Table S6). 



Results 

Molecular phylogenetic analysis of S. kowalevskii Fox 
proteins 

Comparative studies among bilaterian lineages have 
previously identified 24 Fox families in bilaterians (FoxA 
to FoxS), including the newly identified FoxAB family 
[19-21,103]. Two of these families are vertebrate-specific 
(FoxR and FoxS) [29,104], leaving 22 Fox families that 
were present ancestrally in bilaterians. The S. kowalevskii 
genome contains a single copy of all bilaterian Fox family 
members with two notable exceptions; we failed to iden- 
tify a representative of FoxH, which was likely secondarily 
lost in hemichordates and echinoderms [20] since it is 
present in the mollusc Lottia gigantea [33]. Further, the 
FoxQ2 family is represented by three paralogs in S. 
kowalevskii (Figure 1). All S. kowalevskii Fox proteins 
group into their corresponding families with high boot- 
strap and posterior probability values, most closely related 
to the sea urchin or cephalochordate orthologs. Two of 
the sea urchin Fox proteins used in this study are long 
branching (SpFoxX and SpFoxY), and in our analysis they 
could not be assigned to a specific Fox family, supporting 
the result from Tu et al. [20]. The presence of a member 
of the FoxE family in S. kowalevskii, which could not be 
identified in the sea urchin genome, suggests a loss of 
FoxE somewhere in the echinoderm lineage. The FoxAB 
ortholog identified in S. kowalevskii groups reliably together 
with putative orthologs from sea urchin [20], cephalochor- 
dates [21], and cnidarians [19] thus making us confident 
that we have identified an additional member of this new 
family absent in vertebrates (see Figure la, Additional file 7: 
Figure SI and Additional file 8: Table S7). 

FoxQ2 family evolution 

Many animals have multiple duplications in the FoxQ2 
family that have been considered species-specific dupli- 
cations. Multiple FoxQ2 genes are present in the cnidar- 
ians Nematostella vectensis (NvFox4, NvFoxQ2b, NvFoxQ2c) 
[17,63], Hydra magnipapillata (HmFoxQ2b, HmFoxQ2al, 
HmFoxQ2a2) [96], Clytia hemisphaerica (CheFoxQ2a, 
CheFoxQ2b) [19], the cephalochordate Branchiostoma 
floridae (BfFoxQ2a, BJFoxQ2b, BJFoxQ2c) [21], the mollusc 
Lottia gigantea (foxQ2-l, foxQ2-2, foxQ2-3, foxQ2-4), 
and the hemichordate Saccoglossus kowalevskii (foxQ2-l 
[93], foxQ2-2, foxQ2-3). Other animals like the sea urchin 
Strongylocentrotus purpuratus [20], the ascidian Ciona intes- 
tinalis [59], the fish Danio rerio [105], and the fly Drosophila 
melanogaster [37] seem to have only one FoxQ2 gene. 

An ancestral subdivision of the FoxQ2 family at the 
base of the cnidarians has been suggested by Chevalier 
et al. [19], but a more detailed analysis of the evolution- 
ary history of this family was not possible due to lack of 
bilaterian sequences. To address this question we included 
the newly available data from 5. kowalevskii and other 
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Figure 1 (See legend on next page.) 
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(See figure on previous page.) 

Figure 1 Phylogenetic analysis. (A) Phylogenetic analysis of S. kowalevskii Fox genes: The 5. kowalevskii Fox proteins group into their predicted 
families with high support values. Displayed is the Bayesian tree (standard deviation = 0.0109) with Bayesian posterior probabilities values on top 
of each branch and maximum likelihood values underneath each branch. Stars indicate a different tree topology result from the maximum 
likelihood analysis which lead to no support value at that position. Branches with posterior probabilities below 50% are collapsed. For gene 
accession numbers, gene predictions, and alignment see Additional file 1: Table SI, Additional file 2: Table S2, Additional file 3: Table S3 and 
Additional file 4: Table S4. (B) Phylogenetic analysis of the FoxQ2 family. Phylogenetic analysis of FoxQ2 proteins containing an EH-l-like motif 
(see Additional file 5: Table S5) result in a tree topology supporting a duplication of the FoxQ2 family at the base of the bilaterians. Displayed is 
the Bayesian tree (standard deviation = 0.023) with Bayesian posterior probabilities values on top of each branch and maximum likelihood values 
underneath each branch. Stars indicate different tree topologies which lead to no support value at that position. Branches with posterior 
probabilities below 50% are condensed. Proteins with a C-terminal EH-l-like motif are highlighted in blue. Proteins with a N-terminal EH-l-like motif 
are highlighted in yellow. Proteins with a N-terminal and a C-terminal EH-l-like motif are highlighted in yellow and blue. For gene accession numbers, 
identification of the EH-like motif, and alignment see Additional file 1: Table S1, Additional file 2: Table S2, Additional file 3: Table S3, Additional file 4: 
Table S4, Additional file 5: Table S5 and Additional file 6: Table S6. 



phyla (see Additional file 5: Table S5). In our analysis, we 
found that bilaterian FoxQ2 proteins clustered into two 
well supported groups suggesting the duplication of an an- 
cestral FoxQ2 gene occurred before the origin of the bila- 
terians. However, it is not clear whether this event was in 
stem bilaterians or earlier, before the split of cnidarians 
and bilaterians: two cnidarian sequences cluster within one 
of the bilaterian FoxQ2 groups, whereas the others are 
largely unresolved or demonstrate weak support for group- 
ing into the second bilaterian FoxQ2 group (Figure lb). 
More sequence data from additional groups will be re- 
quired to resolve this ambiguity. 

Since the conserved region of Fox proteins is relatively 
short and shows little sequence variability, we found fur- 
ther support for our results by mapping an additional 
character onto the tree: the position of the EH-I-like 
Groucho binding domain. This domain is found in sev- 
eral Fox families, including the FoxQ2 family [105,106]. 
The EH-I-like motif is either located at the C-terminus 
or at the N-terminus of the FoxQ2 protein [105] outside 
the fork head box. Since its sequence is not included in 
the alignment for our phylogenetic analysis, its position 
can be used as an independent character to analyze the 
evolution of this protein family. 

We identified the eight amino acid long EH-I-like 
Groucho binding motif [106,107] for the S. kowalevskii 
FoxQ2 family by manual sequence alignment and NCBI 
Protein BLAST (see Figure lb and Additional file 5: 
Table S5). We found that one of the two bilaterian 
FoxQ2 groups contains all bilaterian FoxQ2 proteins 
that have the EH-I-like motif at the N-terminus of their 
proteins and the second group contains only bilaterian 
FoxQ2 proteins that have the EH-I-like motif at the 
C-terminus. Protein sequences from bilaterian animals 
with multiple FoxQ2 genes, like Saccoglossus kowalevskii, 
Branchiostoma floridae, and Lottia gigantea, were divided 
up into both groups. 

Our data therefore show that the FoxQ2 family was 
already divided into two distinct groups at the base of 
the bilaterians. Since the data from cnidarians do not 



fully resolve timing of this duplication (see Discussion) 
we currently cannot determine whether the split of the 
FoxQ2 family occurred at the bilaterian base or predates 
the divergence of cnidarians and bilaterians. 

Clustered Fox genes 

Two draft genome assemblies for 5. kowalevskii are cur- 
rently available (Baylor College of Medicin/GenBank: 
ACQM00000000.1 and a HudsonAlpha assembly, Hudso- 
nAlpha Institute for Biotechnology, AL (unpublished data)). 
Using these two drafts, four Fox genes show evidence 
of clustering (Additional file 9: Figure S2); foxC, and 
foxLl are joined on one scaffold and foxQ2-l and foxQ2-3 
are closely linked on a separate scaffold (see Additional 
file 9: Figure S2). In addition, foxF clusters with foxC and 
foxLl depending on the algorithm used (it is linked in the 
BCM assembly but not in the HudsonAlpha assembly). 
Further, we provide evidence of a link of foxQl to the foxF, 
foxC, and foxLl containing scaffold by manual genome 
walking using unassembled trace sequences and by bidir- 
ectional blast of the scaffold ends (see Additional file 10: 
Table S8). However, even though no better match was 
found in the genome, the scaffold ends mostly contain re- 
peats and a final assignment of foxQl and foxF requires 
further characterization. The potential linkage of foxQl, 
foxF, foxC, and foxLl is of particular interest since this 
cluster conservation may be linked to their developmental 
roles in mesoderm development [32,33]. 

Expression analysis 

Saccoglossus kowalevskii is a direct developing enterop- 
neust [74,89-91,108]. Early cleavage is radial forming a 
hollow blastula (Figure 2 (1)), and gastrulation is by in- 
vagination between 16 and 30 hours post fertilization 
(hpf) at 20°C (Figure 2 (2)). Mesoderm forms by entero- 
coely from the archenteron following gastrulation at 
about 36 hpf (Figure 2 (3)). Embryogenesis leads to a tri- 
partite body plan with a prosome/proboscis, a mesosome/ 
collar, and a metasome/trunk (Figure 2 (4)), divided after 
48 hpf by an anterior and posterior collar groove (Figure 2 
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Figure 2 Expression patterns of S. kowalevskii foxA-E. Spatial expression pattern of 5. kowalevskii foxA - foxE. Animals are oriented as indicated 
in cartoons for the corresponding stage if not otherwise specified. For a detailed description of the expression patterns see text. Panels 1 to 5; 
ectoderm = light gray, mesoderm = light blue, endoderm = dark gray, black arrows in panel 4 point at the forming furrows at the boundary 
between the proboscis and collar and collar and trunk respectively. Panels 6, 1 1, 16, 19, 21, and 31 show surface views. Panels 17/18 are lateral 
views. (10) Panel 10 shows a dorsal view, white arrows points at the forming gill pores, the black arrow points at the gap of ectodermal foxA expression 
at the dorsal collar, inlay shows lateral view. (15) Inlay in panel 15 shows ventral view on the mouth opening, white arrow points at the mouth 
opening. (20) White arrow in panel 20 points at the endodermal expression domain of foxB. (30) White star in panel 30 indicates the ectodermal 
expression domain of foxD at the base of the proboscis. The inlay shows a closeup of the posterior gut region. An: animal pole; Veg: vegetal pole; 
L: left; R: right; A: anterior; P: posterior; D: dorsal; V: ventral; ao: apical organ; cb: ciliated band; gp: gill pore. Brightness and contrast of pictures 
were adjusted when appropriate to match overall appearance of the figure. 



(4) black arrows). The mouth opens on the ventral side, 
between the collar and proboscis, into the anterior phar- 
ynx, which leads to the posterior gut. The first gill slit 
forms in the posterior pharynx and perforates through the 



ectoderm in the anterior trunk (Figure 2 (5)), with more 
added sequentially during development. A more detailed 
description of hemichordate development can be found in 
[74,88-91,109,110]. 
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foxA 

foxA expression is first detected at blastula in the vegetal 
plate, which gives rise to the prospective endomesoderm 
(Figure 2 (6)). Expression is persistent in the endoderm 
throughout development (Figure 2 (7-10 and inlay)), but 
excluded dorsolaterally from the regions that give rise to 
the gill pores (white arrows in Figure 2 (10)). foxA is also 
expressed in the ectoderm in the anterior collar groove 
from approximately 48 hpf (Figure 2 (9)). In juveniles, 
this circumferential expression marks the most anterior 
collar region but is excluded from the dorsal midline 
(black arrow in Figure 2 (10)). A partial description of 
the expression of foxA was previously reported [73,93]. 

foxAB 

foxAB expression was not detected at blastula stages 
(Figure 2 (11)). At gastrula foxAB is expressed in a cir- 
cumferential ring in the prospective anterior ectoderm 
(Figure 2 (12)). The ectodermal expression persists into 
later stages and refines into a thin ring in the anterior 
collar groove (Figure 2 (12-15)). The developing mouth 
of the embryo perforates through this ring of expression 
on the ventral side (Figure 2 (15 inlay, white arrow indi- 
cates mouth)). 

foxB 

foxB expression is first detected at the blastula stage in a 
circumferential ring in the most posterior prospective 
ectoderm surrounding the vegetal plate (Figure 2 (16)). 
During gastrulation, foxB expression localizes asymmet- 
rically to the posterior ventral ectoderm, flanking the 
ciliated band on both sides (Figure 2 (17-19, Additional 
file 11: Figure S3)). At 48 hpf foxB is expressed in the 
ventral endoderm in the collar region. It is further 
expressed ventrally in a broad stripe in the trunk ecto- 
derm, anterior to the ciliated band and in two further 
narrower stripes posterior to the ciliated band and in 
the collar (Figure 2 (19)). This expression persists into 
the juvenile stage (Figure 2 (20)). At this stage, the 
ectodermal expression domain anterior to the ciliated 
band is divided into two domains (See Additional file 11: 
Figure S3). 

foxC 

foxC expression is first detected at the blastula stage in 
the vegetal plate (Figure 2 (21)). During gastrulation this 
endomesodermal expression restricts to the tip of the 
archenteron, which is fated to become the anterior 
mesoderm [111]. Circumferential ectodermal expression 
is also detected in the anterior of the embryo at the base 
of the developing proboscis (Figure 2 (22)). Upon com- 
pletion of gastrulation, foxC is further associated with 
mesoderm formation and is localized to two pairs of lateral 
endomesoderm (Figure 2 (23)) that become the coelomic 



pouches of the collar and trunk [111] in a pattern very 
similar to foxF (see below). At this stage, circumferential 
ectodermal expression is detected in sharply delineated 
rings marking the anterior and posterior collar groove 
(Figure 2 (23)). From 48 hpf on, foxC expression per- 
sists in the anterior collar groove with a gap of expres- 
sion at the dorsal midline, and endodermal expression 
is detected in the pharyngeal endoderm and is associated 
with gill pouch formation (Figure 2 (24/25)) (See also 
Additional file 11: Figure S3). 

foxD 

foxD expression begins during gastrulation in an anterior 
ectodermal circumferential ring (Figure 2 (26/27)). The 
circumferential ectodermal expression domain persists 
throughout development and localizes to the posterior 
proboscis at later stages (Figure 2 (27-30)). Beginning at 
36 hpf, additional ectodermal expression is detected in 
scattered cells in the proboscis and continues into later 
stages. In the endomesoderm, foxD is expressed poster- 
iorly during gastrulation and early axis elongation 
(Figure 2 (28/29)). Following enterocoely of the meso- 
derm foxD is expressed in the posterior-most endo- 
derm that gives rise to the hindgut (Figure 2 (30 inlay)), 
and in the ventrolateral mesoderm (Figure 2 (30)). 

foxE 

foxE expression is first detected at early blastula in the 
prospective ectoderm (See Additional file 11: Figure S3). 
This expression refines to a ring around the animal pole 
before gastrulation (Figure 2 (31)). During gastrulation 
this circumferential expression domain persists (Figure 2 
(32)) and later refines to a narrow circumferential ring 
localized at the base of the prosome (Figure 2 (33)). 
From 48 hpf on, foxE expression is detected in the anter- 
ior pharynx; the region that later gives rise to the stomo- 
chord (Figure 2 (34/35)) (See also Additional file 11: 
Figure S3). 

foxF 

No foxF expression was observed at blastula (Figure 3 
(1)). Expression is first detected during gastrulation in 
the anterior endomesoderm that will give rise to the an- 
terior mesoderm (Figure 3 (2)). At 36 hpf, it is expressed 
in the developing posterior mesoderm before entero- 
coely (Figure 3 (3)), and at 48 hpf in the anterior, mid, 
and posterior mesoderm (Figure 3 (4)) [111]. At later 
stages mesodermal expression is detected laterally on 
both sides of the pharynx (Figure 3 (5 inlay, indicated by 
stars)), in small patches around the posterior gut, in the 
most anterior tip of the embryo directly underlying the 
apical organ (Figure 3 (5 white arrow)) (Additional file 11: 
Figure S3), and the heart/kidney complex at the tip of 
the developing stomocord (Figure 3 (5 and inlay, black 
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Figure 3 Expression patterns of S. kowalevskii foxF-Ll. Spatial expression pattern of S. kowalevskii foxF - foxLl. Animals are oriented as indicated 
in cartoons of Figure 2 (1-5) for the corresponding stage if not otherwise specified. For a detailed description of the expression patterns see text Panels 
4, 8, 9, and 17-19 show surface views. (5) White arrow in panel 5 points at endodermal expression domain of foxF at the tip of the proboscis, black 
arrow points at the heart-kidney complex. Inlay shows dorsal view of the pharynx, black arrow points at the heart-kidney complex. White asterisk 
indicate expression in the pharyngeal mesoderm. (9) Panel 9 shows dorsal surface. Inlay shows ventral surface. (10) Arrow in panel 10 points 
at dorsal mesoderm. Inlay shows dorsal view of the pharynx, black arrow points at dorsal mesoderm. (15) Black arrow in panel 15 points at 
posterior endoderm expression of foxl. Inlay shows dorsal view of the forming gill pores. (20) Inlay shows dorsal view of the forming gill 
pores. Black arrow heads point to gill pouch endoderm. (30) Inlay shows dorsal view of the forming gill pores. An: animal pole; Veg: vegetal 
pole; L: left; R: right; A: anterior; P: posterior; D: dorsal; V: ventral. Brightness and contrast of pictures were adjusted when appropriate to 
match overall appearance of the figure. 



arrow)). Mesoderm expressing foxF during juvenile devel- 
opment always lines the endoderm (visceral mesoderm) 
except for the expression at the tip of the proboscis meso- 
derm (See also Additional file 11: Figure S3). 

foxG 

No expression of foxG was detected during the blastula 
stage (Figure 3 (6)). Expression is first detected in the 
posterior ectoderm during gastrulation in a few single 



cells (Figure 3 (7)). From 36 hpf on, it is expressed in a 
circumferential ring in the anterior third of the embryo 
and in scattered cells throughout the ectoderm (Figure 3 
(8)). At 48 hpf, the ring of expression refines into two 
separate rings of cells at the base of the proboscis but 
does not extend to the dorsal midline (Figure 3 (9 dorsal 
view, inlay shows ventral view)). Single cell expression is 
detected throughout the proboscis ectoderm with a 
greater density at the proboscis tip. In juveniles, foxG is 
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detected in the dorsal mesoderm overlaying the pharynx 
and gut (Figure 3 (10, black arrow, inlay shows dorsal 
view)). A partial description of foxG expression was pre- 
viously published [71,73]. 

fox/ 

Localized foxl expression was not detected at early de- 
velopmental stages (Figure 3 (11-13)) and is first de- 
tected at 48 hpf in the most posterior endoderm and in 
an ectodermal circumferential ring of expression anter- 
ior to the ciliated band (Figure 3 (14)). At juvenile 
stages, expression is detected in a small domain in the 
center of the gill pouch endoderm (Figure 3 (15, inlay 
dorsal view)) and in the posterior gut (Figure 3 (15, 
black arrow)) (See also Additional file 11: Figure S3). 

fox/7 

No foxjl expression was detected in blastulae (Figure 3 
(16)). foxjl is expressed from gastrula in an ectodermal 
ring marking the forming ciliated band (Figure 3 (17)). 
The expression in the ciliated band is persistent in all 
later stages examined (Figure 3 (17-20)). In juveniles, 
additional broad ectodermal expression is detected in 
the proboscis and the anterior trunk (Figure 3 (20)) in- 
cluding the apical organ (Figure 3 (20 white arrow head)). 
Endodermal expression is detected in the gill pouches 
(Figure 3 (20 and inlay, black arrow heads)) (see also 
Additional file 11: Figure S3). 

foxK 

Only weak staining of foxK was observed at the blastula 
stage (Figure 3 (21)). An almost ubiquitous expression of 
foxK is detected at the late gastrula stage throughout the 
ectoderm but excluded from the ciliated band (Figure 3 
(22)). This expression persists throughout development 
(Figure 3 (22-24)) until it is restricted to the collar ecto- 
derm in juvenile embryos (Figure 3 (25)). Endoderm and 
mesoderm expression was not examined in early devel- 
opmental stages, and no mesoderm or endoderm stain- 
ing was observed in juveniles. 

foxZ.7 

No foxLl expression was detected in early development 
(Figure 3 (26-28)). foxLl expression is first detected at 
48 hpf in a circumferential ectodermal band at the base of 
the proboscis and in the endoderm of the out-pocketing 
gill pouches and in all subsequent stages examined 
(Figure 3 (29/30 inlay shows dorsal view)). 

foxA/7/4 

Ubiquitous expression of foxNl/4 is observed at blastula 
(Figure 4 (1)). At 24 hpf and 36 hpf foxNl/4 is expressed 
throughout the ectoderm except the ciliated band. Scat- 
tered cells nested within the broad domains of ectodermal 



expression show higher levels of expression compared to 
the general ectodermal staining (Figure 4 (2/3)) (See also 
Additional file 11: Figure S3). At 48 hpf, foxNl/4 expres- 
sion continues in the anterior ectoderm but only faint 
staining was detected in the central and posterior collar 
region (Figure 4 (4)). At juvenile stages, foxNl/4 expres- 
sion is restricted to the proboscis and collar, and to a thin 
row of cells (black arrow) posterior to the ciliated band 
(white arrow) (Figure 4 (5, inlay shows ventral view of 
the trunk tip)). Endoderm and mesoderm expression 
was not examined. (For surface views see Additional 
file 11: Figure S3). 

foxP 

Ubiquitous ectodermal foxP expression is detected at 
blastula (Figure 4 (6)). From gastrulation, expression is 
detected in the entire ectoderm except the ciliated band 
(Figure 4 (7)), and stronger expression is detected in sin- 
gle cells scattered throughout the anterior (Figure 4 (7-9 
white arrows)). From 48 hpf, expression is detected in 
single cells throughout the proboscis ectoderm (Figure 4 

(9) ), this expression persists into the juvenile stage (Figure 4 

(10) ). Endoderm and mesoderm expression was not 
examined. 

foxQ7 

No foxQl expression was detected up to 36 hpf of devel- 
opment (Figure 4 (11-13)). Localized expression is first de- 
tected at 48 hpf in the anterior endoderm at the position 
where the gill pouches are forming (Figure 4 (14/15)). At 
juvenile stage, expression continues in the anterior pharynx 
(Figure 4 (15)). If embryos are stained long expression in 
the overlying ectoderm of the trunk and collar, with the ex- 
clusion of the ciliated band, is also observed (Additional 
file 11: Figure S3). 

foxQ2 

All three FoxQ2 paralogs (foxQ2-l, foxQ2-2, and foxQ2-3) 
share an apical expression pattern, but each exhibits unique 
expression characteristics (Figure 4 (16-30)). foxQ2-l is 
expressed in the animal hemisphere at blastula (Figure 4 
(16)). During development from gastrula to juvenile, ex- 
pression becomes refined to the most apical cells marking 
the site of the ciliated apical organ (Figure 4 (17-20)). A 
partial description of the expression pattern of foxQ2-l 
was reported previously by Darras et al. [93] . foxQ2-2 ex- 
pression also displays an ectodermal apical domain similar 
to foxQ2-l (Figure 4 (21-25)). However, in contrast to 
foxQ2-l,foxQ2-2 shows ubiquitous ectodermal expression 
throughout the embryo if stained for longer period of time 
(Additional file 11: Figure S3). Localized foxQ2-3 expres- 
sion is first detected at the gastrula stage in the apical ter- 
ritory (Figure 4 (27)). At 36 hpf, expression in the apical 
domain extends in a stripe along the dorsal midline of the 
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Figure 4 Expression patterns of S. kowalevskii foxN1/4-Q2-3. Spatial expression pattern of S. kowalevskii foxNl/4 - foxQ2-3. Animals are oriented as 
indicated in cartoons of Figure 2 (1 -5) for the corresponding stage if not differently specified. For a detailed description of the expression patterns see 
text. Panels 4, 6, 29, and 30 show surface views. Panels 14/15 and 23-25 show light stained embryos. For longer stained embryos see Additional file 1 1: 
Figure S3. (5) Black arrow points at expression domain of foxNl/4 in the ventral ectoderm at the posterior tip of the trunk, white arrow points at the 
ciliated band. Inlay shows ventral view of the posterior tip of the trunk, black arrow points at expression domain of foxNl/4. (7-9) White arrows 
point to cells with high levels of foxP expression in the proboscis ectoderm. An: animal pole; Veg: vegetal pole; L: left; R: right; A: anterior; P: posterior; 
D: dorsal; V: ventral. Brightness and contrast of pictures were adjusted when appropriate to match overall appearance of the figure. 



embryo starting from the most apical part of the embryo 
and extending posteriorly, covering approximately two- 
thirds of the embryo (Figure 4 (28)). At later stages, this 
dorsal stripe becomes restricted anteriorly to the dorsal 
proboscis midline (Figure 4 (29/30)). 

foxJ2/3, foxL2, foxM, foxN2/3, and foxO 

For five genes foxJ2/3, foxL2, foxM, foxN2/3, and foxO, reli- 
able localization was not detectable by in situ hybridization. 

For a comprehensive summary of all localized S. 
kowalevskii Fox genes see Figure 5. 



Discussion 

Fox gene complement of Saccoglossus kowalevskii 

The Fox gene family is a large group of transcription fac- 
tors with at least 24 family members. Our characterization 
of the Fox gene complement in S. kowalevskii revealed 21 
of the 22 Fox family members predicted for the ancestral 
bilaterian. This includes a member of the FoxE family, 
which is absent in sea urchins, suggesting a loss of FoxE 
somewhere in the echinoderm lineage. S. kowalevskii only 
lacks a representative of the FoxH family, a gene that is 
absent from the genome of the sea urchin S. purpuratus 
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Figure 5 Expression summary. (A-l) Expression summary of all S. kowalevskii Fox genes with clear localized expression patterns. 
(A-E) Blastula -, Gastrula -, 36 h embryo -, 48 h embryo -, 72 h embryo - surface view. (F-H) Gastrula -, 36 h embryo -, 48 h embryo - cross section. 
(I) Cross section of the gill pore area of a 72-h-old embryo. For details see text. *Potential co-expression is inferred from single gene expression analysis. 
No double in situ hybridization was performed. **The expression of foxF is very dynamic and only a more detailed analysis will be able to show all expression 
domains at any given developmental time point. An: animal, Veg: vegetal, A: anterior, P: posterior, D: dorsal, V: ventral. 



and likely secondarily lost in the ambulacrarian lineage 
(See Results and Figure 1). We further identified three 
FoxQ2 genes in S. kowalevskii. The FoxQ2 family likely 
arose in the lineage leading to the common cnidarian bila- 
terian ancestor since no FoxQ2 homologs have been de- 
scribed in sponges [17,33,66] or ctenophores [112]. From 
our phylogenetic analysis, we have identified a subdivision 
of the FoxQ2 family into two distinct groups that occurred 
at least at the base of the bilaterians, but possibly earlier 
(Figure lb and Results). This interpretation is supported 
by the position of the EH-I binding motif, which is located 
either on the C-terminal or N-terminal side of the bilater- 
ian FoxQ2 proteins. Each of the two bilaterian FoxQ2 
groups is characterized by either the C-or N-terminal 
motif. All cnidarian FoxQ2 genes with one exception are 
characterized by a C-terminal EH-I motif. Two of these 
group with the bilaterian FoxQ2 group characterized by 
the EH-I motif at the N-terminus. It is therefore likely 
that the position of the EH-I motif was ancestrally at 
the C-terminus and was relocated to the N-terminal 
domain in one copy of the bilaterian FoxQ2 paralogs 
at the base of the bilaterians. The presence of an N-terminal 
motif in the Hydra gene HmFoxQ2a, which is long branch- 
ing and could not be placed in our phylogenetic analysis, ap- 
pears to be secondarily derived. 

The question remains whether a duplication in the 
FoxQ2 family occurred at the bilaterian base or before 
the bilaterian/cnidarian split. The latter scenario is the 
most parsimonious based on our analysis, but more data 
from cnidarians are needed to answer this question 



definitively: most available cnidarian sequences cur- 
rently group outside all bilaterian FoxQ2 sequences with 
only a few inside the N-terminal clade making this a diffi- 
cult issue to resolve. 

Conserved expression domains of Fox genes 

By comparing the results of our expression analyses 
(Figures 2, 3, 4, and 5) to the current literature we highlight 
the expression of several Fox genes that may represent an- 
cestral deuterostome or bilaterian developmental roles. 

Apical ectoderm patterning and tissue specification 
(foxQ2) 

FoxQ2 gene expression has been characterized during 
the development of a phylogenetically wide range of 
bilaterian and non-bilaterian animals (For references see 
Result section for FoxQ2 and Table 1). In bilaterians 
FoxQ2 genes are commonly expressed at the animal 
pole early in development and quickly restricted to a 
narrow region at the most apical region of the embryo 
[20,21,37,61]. Since this pattern is observed in many 
bilaterians it is proposed that FoxQ2 has an evolution- 
arily conserved role in apical ectoderm patterning [45,61]. 
Functional studies in bilaterians further imply that the re- 
striction of FoxQ2 to the apical tip of the ectoderm is me- 
diated by Wnt/|3-catenin signaling [45,93,113] and imply 
that this regulation is also evolutionarily conserved. 

Recent functional studies outside bilaterians also dem- 
onstrate a regulatory interaction between Wnt/p-catenin 
signaling and FoxQ2 expression [19,65]. In the cnidarian 
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Table 1 Literature summary 



Gene name 



Observed expression in S. kowalevskii 



Related gene expression domain in other 
organisms 



foxA 



Late blastula: 

Vegetal plate. 
48 hpf: 

Anterior collar groove ectoderm; entire endoderm. 
Juvenile stage: 

Entire endoderm with the exception of the 
dorsolateral outpocketing gill pores; most anterior 
collar ectoderm with a gap in expression on the 
dorsal midline 



Protostomes: 

Endoderm specification [103,114-116]. 
Hemichordates: 

Vegetal plate, endoderm/foregut [1 17], 
Echinoderms: 

Presumptive ventral ectoderm [1 18]; Endomesoderm 
specification [44,46]. 

Vertebrates: 

Gastrulation, endoderm patterning, notochord 
formation [18,119]. 

Urochordates: 

Gastrulation, axis formation [1 20]. 
Cnidarians: 

Preendodermal plate, pharynx [62,121]. 



foxAB 



24 hpf: 

Circumferential ectodermal ring which localizes to 
the anterior collar groove during later development. 

Juvenile stage: 

The developing mouth of the embryo breaks through 
at the ventral side of the embryo inside the foxAB 
expression domain 



Protostomes: 

In the bryozoans larval ciliated cleft and apical 
ectodermal territory [103]. 

FoxAB orthologues are further identified in 
echinoderms [20], cephalochordates [21], 
cnidarians [19] but expression is unknown. 



foxB 



Late blastula: 

Ring pattern surrounding the vegetal plate. 
48 hpf: 

Ventral ectoderm anterior and posterior to the ciliated 
band. 

Juvenile stage: 

Ventral ectoderm and ventral pharyngeal endoderm. 



Sea urchins: 

hgressing primary mesenchyme cells; asymmetrically 
expressed along the oral-aboral axis in the oral 
ectoderm and endoderm [122-124]. 

Cnidarians: 

Around the site of gastrulation, larval endoderm 
region, in the statocyst, gonad and tentacle bulb 
of the medusa [19,63]. 

Chordates: 

Dorsal ectoderm, mesoderm, brain: Xenopus [125-127], 
ascidians [59], B. floridae [128], 

Protostomes: 

Ventral ectoderm [40,129]; ventral nervous system [35]. 



foxC 



Late blastula: 

Vegetal plate. 

24 hpf: 

Anterior endomesoderm. 
48 hpf: 

Circumferential ectodermal expression in the anterior 
and posterior collar groove; endomesodermally at the 
positions of mesoderm formation. 

Juvenile stage: 

Endodermal expression in the pharyngeal endoderm 
associated with gill pouch formation. 



Chordates: 

Pharyngeal mesoderm patterning [31,32,34,130-132]. 
Pharyngeal endoderm [34]. 

Vertebrates: 

Ventral and lateral mesoderm, lateral border of 
neuroectoderm, eye, pronephros [133,134] (Xenopus, 
reviewed in [122]) 

Cnidarians: 

Presumptive endoderm and mesenteries [63], 
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Table 1 Literature summary (Continued) 



foxD 



24 hpf: 

Anterior ectodermal circumferential ring. 
36 hpf: 

Anterior ectodermal circumferential ring. 
Posterior endomesoderm. 
48 hpf: 

Anterior ectodermal circumferential ring. 
Ectodermal single cells throughout the proboscis. 
Juvenile stage: 

Posterior-most endoderm forming the hindgut; 
ventrolateral mesoderm. 



Chordates: 

Notochord [52,57] 

Diencephalon: reviewed in [61] 

Neural crest [52,76,135-142] (Reviewed in [52,122]). 
Maintenance of dorsolateral mesoderm (Xenopus) 
[130,143] (reviewed in [122]). Hindgut [135,144], 

Urochordates: 

Dorsal anterior ectoderm [57], 
Cephalochordates: 

Anterior neural plate, the anterior somites, the 
neural tube, and later in the cerebral vesicle, 
hindgut [52,122,145-149], 

Sea urchins: 

Hindgut [20], 
Protostomes: 

Dorso-ventral circumferential cell migration and 
axon projection; ventral mediolateral muscles, 
intestinal precursor cells (C elegans) [150-153], 
Ventral nervous system {Drosophila) [19], 

Cnidarians: 

Aboral third of the embryo, tentacle buds [63]. 
Interaction with BMP/TGF-beta pathway: 

[36,40,55,1 30,143,1 50,1 5 1 ,1 54,1 55]. 



foxE 


12 hpf: 


Vertebrates: 




Ectodermal ring around the animal pole. 


Thyroid (Endoderm) [156] 




48 hpf: 


Urochordates: 




Ectodermal circumferential ring localized at the base 
of the prosome (48 hpf). 


Endostyle (Endoderm) [58] Reviewed in [157]. 




Juvenile stage: 


Cephalochordates: 




Anterior-dorsal pharynx endoderm including the 
stomochord. 


Club shaped gland (Endoderm) [51,158]. 


foxF 


24 hpf: 


Chordates: 




Anterior endomesoderm. 


Mesoderm patterning [19,31,32,34,159-161], 




36 hpf: 


Gill slit formation in chordates: reviewed in [34]. 




Developing lateral and posterior mesoderm. 


Neural plate border, cephalic neural crest, 
pronephros: Xenopus, reviewed in [122], 




48 hpf: 


Protostomes: 




Posterior, central, and anterior mesoderm. 


Mesoderm [36], 




Juvenile stage: 






Mesoderm surrounding the pharynx; mesoderm 
around the posterior gut; a mesodermal spot 
underneath the site of apical organ formation; 
heart-kidney complex; the pharyngeal mesoderm 
with accumulation of foxF expressing cells anterior 
and posterior to the forming gill pores. (Expression 
is absent at the position where the evaginating 
gill pore endoderm connects to the ectoderm.) 
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Table 1 Literature summary (Continued) 


foxG 


24 hpf: 


Vertebrates: 




Few single cells in the ectoderm. 


Telencephalon (reviewed [50]) [162-164]. 




36 hpf: 


Cephalochordates: 




Strong ectodermal circumferential ring in the 
anterior third of the embryo. 


Scattered cells surrounding the cerebral vesicle and 
inside the cerebral vesicle [50] 




48 hpf: 






Two ectodermal rings with a gap of expression 
on the dorsal midline; single cell expression 
throughout the proboscis with a density of 
single cells at the proboscis tip. 






Inx/f^nilp ct^np* 

JUVCIIIIC 31CIUC 






Additionally to 48 hpf expression: dorsal 
mesoderm overlaying the pharynx and gut. 




foxH 


Not present in the 5. kowolevskii genome 




foxl 


48 hpf: 


Vertebrates: 




Most posterior endoderm; weak ectodermal 
circumferential ring of anterior to the ciliated 
band. 

Juvenile stage: 


Craniofacial development [122,127,165]. 
Sea urchins: 

Larval hindgut with high expression levels on the 
aboral side [20]. 




Small domain in the center of the outpocketing 
gill pouch endoderm; posterior gut. 




foxJI 


24 hpf: 


Vertebrates: 




Ectodermal in the ciliated band domain. 
Juvenile stage: 

Ectodermal in the ciliated band, anterior proboscis 
ectoderm including the apical organ; gill pores; 
posterior collar. 


Master regulator in the formation of motile cilia 
[60,166-172], 

Mediates left-right asymmetry [166,173-177]. 
Echinoderms: 

Oral side of the apical plate [20]; Larval ciliary 
band [48]. 

Protostomes: 

Ampullary cells, crescent cells, and prototroch 
(Platynereis) [1 1 3] 

Cnidarians: 

[17,178] 
Yeast: [1 78] 

(Hcm1 p) is involved in spindle pole body 
formation [179], 

FoxJI orthologues are further identified in 
Choanoflagellates [180], sponges [66], other 
deuterostomes and protostomes [33], but 
expression and function is not known. 


foxJ2/3 


Not determined. 




foxK 


>12 hpf: 


Vertebrates: 




Ubiquitous throughout the ectoderm with the 
exception of the ciliated band. 


Dorsal midline, lateral cephalic neural crest, brain, eye, 
lateral muscle precursors {Xenopus) [122,168]. 




Juvenile stage: 






Collar ectoderm. 




foxL1 


48 hpf: 


Chordates: 




Circumferential ectodermal band in the anterior 
collar groove and in the endoderm of the outpocketing 
gill pores. This expression persists until the one gill slit 
stage. 


Pharyngeal mesoderm patterning [31,32,34]. 
Gill slit endoderm (Scyliorhinus canicuia) [34]. 
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Table 1 Literature summary (Continued) 



foxL2 


Not determined. 




foxM 


Not determined. 




foxN1/4 


1 2 hpf: 


Mammals: 




Ubiquitous expression in the ectoderm. 
24 hpf: 


foxNl is essential for proper immune response in mice 

r 1 C5 1 "l DrHA/nstrpsm tsrnpt nf thp Wnt-nsthw^v F1821 

[_ 1 (J 1 J, 1 'UVV 1 1 J LI CO 1 1 1 LdlLJCL LjI UIC V V 1 1 L kJQ L 1 1 V VQ V [ 1 UZJ , 




Ubiquitous ectodermal expression with the 
exception of the ciliated band. 

48 hpf: 

Entire ectoderm with the exception of the central 
and posterior collar region and the ciliated band. 
Nested inside these expression domains are single 
cells that show a high levels of expression. 


foxN4 is involved in specifying amacrine and horizontal 

eel s in the retina and is upstream of the bHLH gene Math3, 

NeuroD1,and Proxl [113,183,184] (reviewed in [185]). 

It is necessary for the development of V2a and V2b 
interneurons in the spinal cord using lateral inhibition 
via the Delta-Notch pathway by activating the transcription 
of Delta4 and the bHLH gene Mash-1 [186,187], 




Juvenile stage: 






Proboscis and collar region; thin row of cells 
posterior to the ciliated band. 




foxN2/3 


Not determined 




foxO 


Not determined 




foxP 


24 hpf: 


Vertebrates: 




Entire ectoderm with the exception of the ciliated 
band. 

>36 hpf: 

High levels of expression in single cells in the 
anterior ectoderm of the embryo; weak 
endodermal expression. 


Basal ganglia, cerebral cortex, cerebellum, and thalamus, 
hippocampus [127,188-190]. In Medaka foxPl expression 
indicates a role during striatum projection neuron 
development [191] and basal ganglia development of 
the developing and central nervous system [192]. 

Mutations in human foxP2 gene lead to severe 
language disorders [193-195] (reviewed in [127,196]). 

Urochordate: 

Developing brain [59]. 

FrhinnHnrmc 

111 IUUCI 1 1 ID. 

Fore- and mid-gut of the larva [20]. 
Protostomes: 

In Drosophila the two splice variants of foxP (fd85Ea 

and fd85Eb) are expressed in the developing CNS [37,197]. 

A FoxP ortholog is also found in other ecdysozoans, 
cnidarians, and sponges [33,66], but expression 
patterns are not yet reported. 


foxQ1 


48 hpf: 


Vertebrates: 




Anterior endoderm at the position where the gill 
pores are forming; faint circumferential ectodermal 
ring at the position where the gill pores are forming. 

Juvenile stage: 

Anterior pharynx. 


Prospective pharynx, pharyngeal pouches [60], 

Cephalochordates: 

Endostyle [54], 

Urochordates: 

Pharyngeal gills, endostyle [58]. 
Protostomes: 

Pharyngeal endoderm [32]. 


foxQ2-1 


1 2 hpf: 


Deuterostomes: 




Animal hemisphere. 


Apical ectoderm [20,21,53,198]. 




>24 hpf: 


In protostomes: 




During the development from gastrula to juvenile, 
the ectodermal expression becomes refined to cells 
forming the apical organ. 


Anterior tip of the embryo at the blastoderm stage, 
pharyngeal structures and the brain hemispheres 
{Drosophila) [37]. Apical ectoderm (Platynereis) [1 1 3] 
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Table 1 Literature summary (Continued) 



Cnidarians: 

Apical tuft [19,178], 

Linke to Wnt/p catenin signaling: 

[1 9,45,64,73]. 



foxQ2-2 



foxQ2-2 expression resembles the expression 
of foxQ2-1 with additional ubiquitous expression 
throughout the embryo ectoderm. 



foxQ2-3 



24 hpf: 

Apical ectodermal territory. 
36 hpf: 

Ectodermal stripe along the dorsal axis of the 
embryo starting from the most apical part of 
the embryo and extending posterior, covering 
approximately two-thirds of the embryo. 

Juvenile stage: 

Ectodermal stripe along the dorsal proboscis 
midline. 



Clytia hemisphaerica, CheFoxQ2a is expressed opposite 
the site of gastrulation, which is marked by nuclear 
|3-catenin and the expression of multiple Wnt genes 
[19,64]. Knock down of CheWnt3a prevents the restric- 
tion of CheFoxQ2a to the most aboral ectodermal tip of 
the embryo [64,65] demonstrating a regulatory inter- 
action of CheFoxQ2a with the canonical Wnt signaling 
pathway during axial patterning and supporting the hy- 
pothesis of an evolutionarily conserved interaction of 
FoxQ2 and Wnt/|3-catenin signaling predating the split 
between cnidarians and bilaterians. 

In addition to the role of FoxQ2 in apical ectoderm 
patterning it is also proposed to be involved in apical organ 
formation, a neural rich structure located at the apical tip 
of many bilaterian and non-bilaterian larvae [199] generally 
comprised of sensory cells, neurons, and long motile cilia 
that form the apical tuft. Even though a homology of 
the apical organ among metazoans is favored [113,178, 
200-203] there remain dissenting views [61,204,205]. 

In both bilaterian and cnidarian larvae with an apical 
organ/tuft, FoxQ2 expression coincides with the position 
of this structure (see Figure 6) [19,178]. Functional stud- 
ies specifically investigating the role of FoxQ2 in apical 
organ formation have been demonstrated in sea urchins 
and cnidarians [178,206]. Loss of function in sea urchins 
compromises the development of the long apical tuft 
cilia [206]. In the cnidarian N. vectensis expression of 
one out of four FoxQ2 genes has been reported [178]. 
At the planula stage it is expressed around the apical 
organ/apical tuft. Knock down experiments show that 
NvFoxQ2a is involved in the determination of the size of 
the apical organ/apical tuft [178]. 

Expression data and functional studies further suggest 
an evolutionarily conserved core regulatory network for 



apical organ formation including an apical Six3 domain, 
with FoxQ2 and Foxjl expressed in the apical organ do- 
main in a Wnt negative territory [113,178]. 

In S. kowalevskii expression of all three FoxQ2 genes is 
localized at the apical tip of the ectoderm during develop- 
ment. foxQ2-l, which shows the most resemblance in ex- 
pression to other reported FoxQ2 expression patterns, 
begins broadly in the animal hemisphere at early blastula 
stage. At later developmental stages, expression becomes 
restricted to the apical tip at the site of the apical organ 
(Figure 4 (16-20)) similar to other bilaterian groups. foxQ2- 
1 is co-expressed with the motile cilia marker foxjl in the 
apical organ (see below and Result section for foxjl), but 
unlike other groups foxQ2-l is co-expressed with six3 [71], 
and foxjl has other broad expression domains not tightly 
localized to the apical organ. Experimental manipulations 
of Wnt/[3-catenin signaling in S. kowalevskii provide sup- 
port for a Wnt/|3-catenin dependent localization of foxQ2-l 
to the apical ectoderm [93], but a role of foxQ2-l in apical 
tuft patterning will still need to be directly tested by ex- 
perimental knockdown. Our analysis of S. kowalevskii 
FoxQ2 genes therefore supports the hypothesis that they 
play an evolutionarily conserved role in patterning an an- 
terior ectodermal territory in bilaterians and that an- 
terior restriction of FoxQ2 is mediated by Wnt/[3-catenin 
signaling. 

Formation of motile cilia {foxJI) 

In contrast to primary cilia which have a 9 + 0 arrange- 
ment of their microtubules and only sometimes possess 
dynein arms (Nodal cilia), motile cilia show a 9 + 2 ar- 
rangement of their microtubules, are generally longer 
than primary cilia, and possess outer dynein arms medi- 
ating motility (reviewed in [207,208]). In vertebrates, 
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Figure 6 Examples of conserved Fox gene expression domains. (A-E) FoxB expression (blue) in multiple species. The expression of FoxB 
seems to be correlated to the expression of chordin, a BMP inhibitor. The endodermal gut expression domain seems to be unique to 
echinoderms and hemichordates. (F-l) FoxD expression (blue) in the hindgut of different species. The conservation of the expression in the 
hindgut across deuterostomes indicates that it was likely already present in the hindgut of the deuterostome ancestor. (J-0) Ectodermal 
expression domains of FoxQ2 (blue) in multiple species across phyla. Expression across bilateria suggest a conserved role in apical (apical neural) 
patterning. Expression at the aboral side in cnidarians suggests a deep evolutionary origin for this expression in patterning terminal neural 
structures. For literature summary see Table 1 and for additional discussion see Additional file 12. A: anterior, P: posterior, D: dorsal, V: ventral, 
oa: aboral, o: oral. 



motile cilia are mainly found in tissues where fluid 
movement is necessary, like lung epithelia or in the em- 
bryo node where it is involved in mediating left-right 
asymmetry [166,173-176,209]. Expression and functional 
studies in vertebrates imply that Foxjl is a master regula- 
tor of the formation of these motile cilia [60,166-169,210] 
(reviewed in [209]). A conserved role of Foxjl in motile 
cilia formation is supported from expression patterns out- 
side chordates.^wc/i expression has been described during 
sea urchin development where it is expressed in the most 
apical ectoderm marking the apical organ/apical tuft, and 
the ciliated bands [20,48], the sea star Patiria miniata 
where foxjl is expressed throughout the ectoderm at the 
blastula stage and the ciliated bands at the larval stage 



[48], the planarian Schmidtea mediterranea where Smed- 
foxJl-2 and Smed-foxJl-2 are expressed along the AP axis 
in a ventral stripe of presumptive motile ciliated sensory 
cells [210], the annelid Platynereis dumerilii where foxjl is 
expressed in the apical plate and the ciliated bands, and 
the cnidarian Nematostella vectensis where NvFoxJl is 
marking the apical organ territory [178]. 

In S. kowalevskii foxjl is also expressed in the region 
forming motile cilia; the ciliated band (Figure 3 (17-20)) 
and the gill pore endoderm in juveniles (Figure 3 (20)). 
In juveniles it is further expressed broadly in the anterior 
ectoderm which incorporates the territory of the apical 
organ. Our data therefore support the hypothesis that 
Foxjl has an evolutionarily conserved function in motile 
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cilia formation that predates the bilaterian-cnidarian split. 
Further, our data support an evolutionary conserved 
function of Foxjl in cilia formation in the apical organ 
despite the rather broad anterior ectodermal expression in 
S. kowalevskii when compared to the tightly localized ex- 
pression of Foxjl in the developing apical organs of other 
metazoans. However, since the expression of Foxjl is cor- 
related with all motile ciliated cells with a 9 + 2 arrange- 
ment of their microtubules, Foxjl expression provides 
limited insights into the homology of this structure. 

Anterior endoderm and mesoderm patterning {foxQ1, 
foxF, foxC, foxLI) 

The third germ layer of bilaterians, mesoderm, likely 
evolved at the base of the bilaterians and gives rise to 
many essential components of the bilaterian body. The 
emergence and evolution of mesoderm is therefore of 
special interest to understand bilaterian body plan evolu- 
tion. Several Fox genes, namely FoxF, FoxC, and FoxLI, 
are proposed to have evolutionarily conserved functions in 
patterning distinct mesodermal populations [31,32,211]. An- 
cestral linkage of these three Fox genes, along with FoxQl 
in stem bilaterians, has been proposed to be related to their 
conserved developmental functions [31,32,34,159,211]. The 
analysis from Shimeld et al. [32] suggests that FoxC and 
FoxLI play a conserved role in somatic mesoderm forma- 
tion (mesoderm lining of ectoderm), FoxF in visceral 
mesoderm formation (mesoderm lining of endoderm), 
and FoxQl in anterior gut endoderm formation. 

Our genomic analysis of foxQl, foxF, foxC, and foxLI 
in S. kowalevskii shows a possible link of these four 
genes in the 5 .kowalevskii genome (see Additional file 9: 
Figure S2). Expression studies show that foxLI is only 
expressed in the endoderm in S. kowalevskii but not in 
the mesoderm. Expression of FoxLI in the mesoderm of 
protostomes and other deuterostomes suggests this likely 
represents secondary loss in hemichordates. Expression of 
foxQl in the foregut of S. kowalevskii is consistent with 
the proposed conserved role in anterior endoderm forma- 
tion. Early expression of foxF in the forming mesoderm 
suggests a conserved role in mesoderm patterning, but the 
division of early mesoderm into somatic and visceral terri- 
tories remains to be characterized in 5. kowalevskii. How- 
ever, in later stages (juveniles) foxF marks the mesoderm 
surrounding the gut (Figure 3 (1-5 and Additional file 11: 
Figure S3)) in both the collar and trunk, consistent with a 
role in visceral mesoderm patterning. foxC is expressed in 
all the mesodermal compartments during initial specifica- 
tion and before out-pocketing of coeloms. However, our 
analysis did not detect later expression in the differenti- 
ated mesoderm making it difficult to speculate whether 
this gene plays a conserved role specifically in somatic 
mesoderm development. 



In conclusion, our data provide support for an ances- 
tral bilaterian chromosomal linkage of foxQl-foxF-foxC- 
foxLl, and an evolutionarily conserved role of FoxF in 
visceral mesoderm patterning and FoxQl in foregut pat- 
terning. However, whether the early expression of foxC 
in mesoderm patterning is related to somatic mesoderm 
formation will require further investigation. 

Gill slit patterning {foxC, foxl, foxJI, foxLI) 

Hypotheses of homology of deuterostome pharyngeal 
gill slits have a long history in comparative studies 
[74]. Morphological and molecular studies in hemichor- 
dates support the homology of pharyngeal gills between 
ambulacrarians and chordates [75-79,212,213]. In this study 
we describe the expression of several Fox genes that can 
further contribute to this discussion. In S. kowalevskii foxC, 
foxl, foxjl, and foxLI are expressed in the endoderm of the 
first gill pouch (Figure 2 (24/25), 3 (15), 3 (29/30)). Recent 
molecular data revealed that FoxC, and FoxLI play con- 
served roles in gill slit formation in chordates, with 
conserved expression in the forming gill slit mesoderm 
(reviewed in Wotton et al. [34]). However, endodermal 
expression of foxC and foxLI during gill formation has 
only been described in the dogfish Scyliorhinus canicula 
[34]. Endodermal expression of Foxl during pharyngeal 
pouch development has been described in mice (foxi3 
[165]) and zebrafish (foxil [214,215]). In S. kowalevskii 
endodermal expression of foxl is detected during the de- 
velopment of the first gill pouch (Figure 3 (15)), support- 
ing a conserved role for this gene in deuterostome gill 
formation. Our data suggest that endodermal expression 
of FoxC, FoxLI, and Foxl was involved in patterning the 
pharyngeal gill endoderm in stem deuterostomes, extend- 
ing the analysis of Gillis et al. [79] and further strengthen- 
ing hypotheses of deuterostome gill slit morphological 
homology by reconstructing ancestral gene regulatory net- 
works involved in early deuterostome pharyngeal endo- 
dermal patterning. 

Ventral endoderm, mesoderm, and ectoderm patterning 
(foxB, foxD) 

Basic anatomical comparisons on the relative organization 
of organ systems across the dorsoventral axis of arthro- 
pods and chordates have resulted in hypotheses suggesting 
the equivalence of the dorsal side of chordates and the 
ventral side of arthropods. This so-called dorsoventral axis 
inversion hypothesis has gained molecular support from 
comparative developmental genetics [216,217] (reviewed 
in [218,219]). In S. kowalevskii bmp2/4 and chordin, which 
are involved in mediating DV patterning in bilaterian 
groups, are expressed in the same relative position as pro- 
tostomes during DV patterning; bmp2/4 dorsally, and 
chordin ventrally, which is inverted relative to their ex- 
pression in chordates. These data suggest a molecular 
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inversion of the DV axis after the split of ambulacrarians 
and chordates [83,219]. In our survey, we revealed two 
Fox genes, foxB and foxD, with differential expression 
along the DV axis, further supporting this observation of 
DV inversion. 

In chordates, FoxB and FoxD are consistendy expressed 
dorsally. In Xenopus laevis foxDl,foxD2, and foxD3 are all 
expressed in the dorsal mesoderm [122] and are necessary 
for dorsolateral mesoderm identity [130,143]. In proto- 
stomes FoxD expression is detected in the ventral nervous 
system of D. melanogaster [35], and the ventral mediolat- 
eral muscles of C. elegans [150,151] where it is known for 
its function in dorso-ventral cell migration and axon 
projection [152,153]. FoxB is expressed in the dorsal 
ectoderm and mesoderm in chordates (frogs [122,125], 
ascidians [59], and cephalochordates [128]). In proto- 
stomes, FoxB expression in C. elegans (lin-31) is localized 
to ventral ectodermal cells [129], and in D. melanogaster, 
the two FoxB orthologs (FD4 and FD5) are expressed in 
the ventral nervous system [35]. 

The localized expression of FoxD and FoxB along the 
bilaterian DV axis raises the possibility of a link to the 
BMP/TGF-beta pathway. Limited comparative functional 
studies confirm a link between FoxD and the BMP/TGF- 
beta pathway [55,154,155,220,221]. In chordates {Xenopus), 
foxDl (XBF-2) is downstream of BMP- antagonists like 
Cerberus, Noggin, and Gremlin, and plays a role in 
maintaining dorsolateral mesoderm during gastrulation 
by downregulating BMP-4 [130,143]. In protostomes 
(C. elegans) FoxD (UNC-130) acts as a transcriptional 
repressor and inhibits the expression of UNC-129, a C. 
elegans TGF-beta ortholog [40,150,151]. However, expres- 
sion data from sea urchins and ascidians are inconsistent 
with a link to BMP signaling [20,57]. In S. kowalevskii, 
foxD and foxB are expressed ventrally opposite the side of 
BMP expression [83] (Figure 2 (17-20), 2 (30)) (Figure 5) 
supporting a link to BMP signaling and DV patterning, 
but this will need to be functionally validated. 

Through gut evolution (foxAB, foxD, fox/) 

In our study we found three Fox genes that are expressed 
either in the mouth or hindgut of the embryo. foxAB, is 
expressed in a circumferential ectodermal ring in the an- 
terior collar groove in S. kowalevskii at the position where 
the mouth forms (Figure 2 (12-15)). S. kowalevskii foxD 
and foxl are expressed in the hindgut (Figure 2 (30), 3 
(15)). FoxD is also expressed in the hindgut of several 
other deuterostome species, including the sea urchin 
S. purpuratus, the cephalochordate B. floridae, the 
frog X. laevis, and the fish O. latipes, and D. rerio 
[20,52,135,144,192] (illustrated in Figure 6 F-I). In the 
protostome groups examined, D. melanogaster and C. ele- 
gans, there is no support for a conserved role of FoxD 
in gut pattering [35,150]. Hindgut Foxl expression similar 



to S. kowalevskii has so far only been described in sea ur- 
chins, where foxl is expressed in the larval hindgut with 
strong expression on the aboral side [20] . 

A FoxAB ortholog was also identified in echinoderms 
[20], cnidarians [19], cephalochordates [21], and bryo- 
zoans [103]. Expression however, is only known from 
bryozoans where it is expressed only transiently in larval 
structures in the ciliated cleft and abapical ectodermal 
territory [103]. Further expression analyses will be re- 
quired before evolutionary hypotheses of its role in 
mouth formation can be tested. 

From these data we conclude that FoxD likely plays a 
conserved role in deuterostome hindgut patterning, Foxl 
was likely co-opted into hindgut patterning during ambu- 
lacrarian evolution, and a broader bilaterian role of FoxAB 
in mouth patterning will require additional data from 
other phyla. 

For a literature summary, see Table 1 and for additional 
discussion on expression patterns see Additional file 12. 

Conclusions 

Analyzing the expression patterns of Fox genes in the 
hemichordate Saccoglossus kowalevskii builds on com- 
parative data from echinoderms and chordates and helps 
to reconstruct the evolutionary history of the develop- 
mental roles of this important family of transcription 
factors during deuterostome evolution. Further compar- 
ing these data to available studies from all metazoans 
helps us to construct more robust hypotheses about the 
role of Fox genes as components of evolutionarily con- 
served gene regulatory networks by distinguishing them 
from lineage specific co-option. 

Our sequence analysis demonstrates that all 23 Fox 
genes of S. kowalevskii fall into their respective families. 
It further refines our understanding of the evolution and 
diversification of the FoxQ2 gene family revealing a split 
of this family deep in metazoan evolution. We provide 
evidence for a clustered arrangement of foxQlfoxF- 
FoxC-foxLl in the S. kowalevskii genome, which has 
been proposed to be an ancestral feature of bilaterians. 

From our expression analyses we propose several evo- 
lutionarily conserved expression domains. In multiple 
cases these gene expressions support hypotheses of ana- 
tomical homology between phyla; conserved expression 
of foxC, foxl, foxjl, and foxLl during gill slit formation 
provides additional molecular support for the presence of 
pharyngeal gills in the common deuterostome ancestor; 
mesodermal foxF expression supports an evolutionarily 
conserved role for FoxF in visceral mesoderm patterning; 
expression of foxQl supports an evolutionarily conserved 
role of FoxQl in pharyngeal endoderm patterning; foxjl 
expression supports a conserved role of Foxjl in forming 
motile cilia throughout metazoans; and a conserved role 
of Foxl and FoxD in hindgut patterning in ambulacrarians 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



Page 20 of 25 



provide evidence of a link of foxQl to the foxF, foxC, and foxLl containing 
scaffold by manual genome walking using unassembled trace sequences 
and by bidirectional blast of the scaffold ends (see Additional file 10: 
Table S8). However, even though no better match was found in the 
genome, the scaffold ends mostly contain repeats and a final assignment 
of foxQl and foxF requires further characterization. The S. kowalevskii foxQ2- I 
and foxQ2-3 cluster indicates a species-specific tandem duplication event. 
Red arrows indicate orientation of the genes, black arrows indicate the 
continuation of a scaffold, and distances are given in kilobase pairs 
underneath each cluster. Black line connecting foxQl and the foxF-foxC-foxL 
cluster indicates area of manual genome walking. 

Additional file 10: Table S8. Bridging contigs for foxQl and foxLI-foxF 

contig. 

Additional file 11: Figure S3. Additional stages and views of Fox gene 
expression patterns. 

Additional file 12: Discussing various Fox gene expression patterns 
and their potential evolutionary relevance. 



and deuterostomes is supported by their expression 
pattern in S. kowalevskii, respectively. In other cases, 
expression in S. kowalevskii supports conserved inter- 
actions with signaling pathways such as FoxQ2 and 
Wnt/(3-catenin signaling, and FoxB and FoxD with the 
BMP pathway. 

Further studies of other phyla, particularly the lopho- 
trochozoans, acoels, and cnidarians, are now required to 
broaden the phylogenetic scope of these comparisons. 
Functional studies are further required to confirm the 
proposed interactions of Fox genes with signaling path- 
ways to further elucidate the evolution of this transcrip- 
tion factor family and its roles in embryonic patterning. 

Additional files 
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BCM assembly but not in the HudsonAlpha assembly). Further, we 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



7. Clark KL, Halay ED, Lai E, Burley SK: Co-crystal structure of the HNF-3/fork 
head DNA-recognition motif resembles histone H5. Nature 1993, 
364:412-420. 

8. Li C, Tucker PW: DNA-binding properties and secondary structural model 
of the hepatocyte nuclear factor 3/fork head domain. Proc Natl Acad Sci 
USA 1993, 90:11583-11587. 

9. Kaufmann E, Hoch M, Jackie H: The interaction of DNA with the DNA- 
binding domain encoded by the Drosophila gene fork head. Eur J 
Biochem 1994, 223:329-337. 

10. Kaufmann E, Knochel W: Five years on the wings of fork head. Mech Dev 
1996, 57:3-20. 

1 1. Carlsson P, Mahlapuu M: Forkhead transcription factors: key players in 
development and metabolism. Dev Biol 2002, 250:1-23. 

1 2. Stroud JC, Wu Y, Bates DL, Han A, Nowick K, Paabo S, Tong H, Chen L: 
Structure of the forkhead domain of FOXP2 bound to DNA. Structure 
2006, 14:159-166. 

13. Jurgens G, Wieschaus E, Nusslein-Volhard C, Kluding H: Mutations affecting 
the pattern of the larval cuticle in Drosophila melanogaster. II Zygotic 
loci on the third chromosome Roux's Arch. Dev Biol 1 984, 1 93:283-295. 

14. Costa RH, Grayson DR, Darnell JE: Multiple hepatocyte-enriched nuclear 
factors function in the regulation of transthyretin and alpha 1-antitrypsin 
genes. Mol Cell Biol 1989, 9:1415-1425. 

1 5. Lai E, Prezioso VR, Smith E, Litvin 0, Costa RH, Darnell JE: HNF-3A, a 
hepatocyte-enriched transcription factor of novel structure is regulated 
transcriptionally. Genes Dev 1990, 4:1427-1436. 

16. Weigel D, Jackie H: The fork head domain: a novel DNA binding motif of 
eukaryotic transcription factors? Cell 1990, 63:455-456. 

17. Larroux C, Luke GN, Koopman P, Rokhsar DS, Shimeld SM, Degnan BM: 
Genesis and expansion of metazoan transcription factor gene classes. 
Mol Biol Evol 2008, 25:980-996. 

18. Kaestner KH, Knochel W, Martinez DE: Unified nomenclature for the 
winged helix/forkhead transcription factors. Genes Dev 2000, 14:142-146. 

1 9. Chevalier S, Martin A, Leclere L, Amiel A, Houliston E: Polarised expression 
of FoxB and FoxQ2 genes during development of the hydrozoan Clytia 
hemisphaerica. Dev Genes Evol 2006, 216:709-720. 

20. Tu Q, Brown CT, Davidson EH, Oliveri P: Sea urchin Forkhead gene family: 
phylogeny and embryonic expression. Dev Biol 2006, 300:49-62. 

21. Yu JK, Mazet F, Chen YT, Huang SW, Jung KC, Shimeld SM: The Fox genes 
of Branchiostoma floridae. Dev Genes Evol 2008, 218:629-638. 

22. Frank S, Zoll B: Mouse HNF-3/fork head homolog-1-like gene: structure, 
chromosomal location, and expression in adult and embryonic kidney. 
DNA Cell Biol 1 998, 1 7:679-688. 

23. Bieller A, Pasche B, Frank S, Glaser B, Kunz J, Witt K, Zoll B: Isolation and 
characterization of the human forkhead gene FOXQ1. DNA Cell Biol 2001, 
20:555-561. 

24. Hong HK, Noveroske JK, Headon DJ, Liu T, Sy MS, Justice MJ, Chakravarti A: 
The winged helix/forkhead transcription factor Foxql regulates 
differentiation of hair in satin mice. Genesis 2001, 29:163-171. 

25. Shu W, Yang H, Zhang L, Lu MM, Morrisey EE: Characterization of a new 
subfamily of winged-helix/forkhead (Fox) genes that are expressed in 
the lung and act as transcriptional repressors. J Biol Chem 2001, 
276:27488-27497. 

26. Schubert LA, Jeffery E, Zhang Y, Ramsdell F, Ziegler SF: Scurfin (FOXP3) acts 
as a repressor of transcription and regulates T cell activation. J Biol Chem 
2001, 276:37672-37679. 

27. Katoh M, Katoh M: Identification and characterization of human FOXN5 
and rat Foxn5 genes in silico. Int J Oncol 2004, 24:1339-1344. 

28. Katoh M, Katoh M: Identification and characterization of human 
FOXN6, mouse Foxn6, and rat Foxn6 genes in silico. Int J Oncol 2004, 
25:219-223. 

29. Heglind M, Cederberg A, Aquino J, Lucas G, Ernfors P, Enerback S: Lack of 
the central nervous system- and neural crest-expressed forkhead gene Foxsl 
affects motor function and body weight. Mol Cell Biol 2005, 25:561 6-5625. 

30. Mazet F, Yu JK, Liberies DA, Holland LZ, Shimeld SM: Phylogenetic 
relationships of the Fox (Forkhead) gene family in the Bilateria. Gene 
2003, 316:79-89. 

31. Mazet F, Amemiya CT, Shimeld SM: An ancient Fox gene cluster in 
bilaterian animals. Curr Biol 2006, 16:R314-R316. 

32. Shimeld SM, Boyle MJ, Brunet T, Luke GN, Seaver EC: Clustered Fox genes 
in lophotrochozoans and the evolution of the bilaterian Fox gene 
cluster. Dev Biol 2010, 340:234-248. 



Page 21 of 25 



33. Shimeld SM, Degnan B, Luke GN: Evolutionary genomics of the Fox genes: 
Origin of gene families and the ancestry of gene clusters. Genomics 2010, 
95:256-260. 

34. Wotton KR, Mazet F, Shimeld SM: Expression of FoxC, FoxF, FoxLI, and 
FoxQI genes in the dogfish Scyliorhinus canicula defines ancient and 
derived roles for Fox genes in vertebrate development. Dev Dyn 2008, 

237:1590-1603. 

35. Hacker U, Grossniklaus U, Gehring WJ, Jackie H: Developmental^ regulated 
Drosophila gene family encoding the fork head domain. Proc Natl Acad 
Sci USA] 992, 89:8754-8758. 

36. Perez Sanchez C, Casas-Tinto S, Sanchez L, Rey-Campos J, Granadino B: Dm 
Fox F, a novel Drosophila fork head factor expressed in visceral mesoderm. 
Mech Dev 2002, 111:163-166. 

37. Lee HH, Frasch M: Survey of forkhead domain encoding genes in the 
Drosophila genome: Classification and embryonic expression patterns. 
Dev Dyn 2004, 229:357-366. 

38. Bulow MH, Aebersold R, Pankratz MJ, Junger MA: The Drosophila FoxA 
ortholog Fork head regulates growth and gene expression downstream 
of Target of rapamycin. PLoS One 2010, 5:el 51 71 . 

39. Kalb JM, Lau KK, Goszczynski B, Fukushige T, Moons D, Okkema PG, 
McGhee JD: pha-4 is Ce-fkh-1, a fork head/HNF-3alpha, beta, gamma 
homolog that functions in organogenesis of the C. elegans pharynx. 
Development 1 998, 125:2171-2180. 

40. Hope IA, Mounsey A, Bauer P, Aslam S: The forkhead gene family of 
Caenorhabditis elegans. Gene 2003, 304:43-55. 

41 . Gao X, Wang Z, Martin J, Abubucker S, Zhang X, Mitreva M, Hawdon JM: 
Identification of hookworm DAF-16/FOXO response elements and direct 
gene targets. PLoS One 2010, 5:el 2289. 

42. Yen K, Narasimhan SD, Tissenbaum HA: DAF-16/Forkhead box O 
transcription factor: many paths to a single Fork (head) in the road. 
Antioxid Redox Signal 201 1 , 1 4:623-634. 

43. Hinman VF, Nguyen AT, Cameron RA, Davidson EH: Developmental gene 
regulatory network architecture across 500 million years of echinoderm 
evolution. Proc Natl Acad Sci USA 2003, 100:13356-13361. 

44. Oliveri P, Walton KD, Davidson EH, McClay DR: Repression of mesodermal 
fate by foxa, a key endoderm regulator of the sea urchin embryo. 
Development 2006, 1 33:41 73-41 81 . 

45. Yaguchi S, Yaguchi J, Angerer RC, Angerer LM: A Wnt-FoxQ2-nodal pathway 
links primary and secondary axis specification in sea urchin embryos. Dev 
Cell 2008, 14:97-107. 

46. Hinman VF, Yankura KA, McCauley BS: Evolution of gene regulatory 
network architectures: examples of subcircuit conservation and 
plasticity between classes of echinoderms. Biochim Biophys Acta 2009, 
1789:326-332. 

47. McCauley BS, Weideman EP, Hinman VF: A conserved gene regulatory 
network subcircuit drives different developmental fates in the vegetal 
pole of highly divergent echinoderm embryos. Dev Biol 2010, 340:200-208. 

48. Yankura KA, Martik ML, Jennings CK, Hinman VF: Uncoupling of complex 
regulatory patterning during evolution of larval development in 
echinoderms. BMC Biol 2010, 8:143. 

49. Shimeld SM: Characterisation of amphioxus HNF-3 genes: conserved 
expression in the notochord and floor plate. Dev Biol 1997, 183:74-85. 

50. Toresson H, Martinez-Barbera JP, Bardsley A, Caubit X, Krauss S: Conservation 
of BF-1 expression in amphioxus and zebrafish suggests evolutionary 
ancestry of anterior cell types that contribute to the vertebrate telencephalon. 
Dev Genes Evol 1998, 208:431-439. 

51 . Yu JK, Holland LZ, Jamrich M, Blitz IL, Holland ND: AmphiFoxE4, an 
amphioxus winged helix/forkhead gene encoding a protein closely 
related to vertebrate thyroid transcription factor-2: expression during 
pharyngeal development. Evol Dev 2002, 4:9-15. 

52. Yu JK, Holland ND, Holland LZ: An amphioxus winged helix/forkhead 
gene, AmphiFoxD: insights into vertebrate neural crest evolution. Dev 
Dyn 2002, 225:289-297. 

53. Yu JK, Holland ND, Holland LZ: AmphiFoxQ2, a novel winged helix/ 
forkhead gene, exclusively marks the anterior end of the amphioxus 
embryo. Dev Genes Evol 2003, 213:102-105. 

54. Mazet F, Luke GN, Shimeld SM: The amphioxus FoxQI gene is expressed 
in the developing endostyle. Gene Expr Patterns 2005, 5:313-315. 

55. Onai T, Yu JK, Blitz IL, Cho KW, Holland LZ: Opposing Nodal/Vgl and BMP 
signals mediate axial patterning in embryos of the basal chordate 
amphioxus. Dev Biol 2010, 344:377-389. 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



Page 22 of 25 



56. Jeffery WR, Ewing N, Machula J, Olsen CL, Swalla BJ: Cytoskeletal actin 
genes function downstream of HNF-3beta in ascidian notochord 
development. Int J Dev Biol 1998, 42:1085-1092. 

57. Imai KS, Satoh N, Satou Y: An essential role of a FoxD gene in notochord 
induction in Ciona embryos. Development 2002, 129:3441-3453. 

58. Ogasawara M, Satou Y: Expression of FoxE and FoxQ genes in the 
endostyle of Ciona intestinalis. Dev Genes Evol 2003, 21 3:41 6-41 9. 

59. Imai KS, Hino K, Yagi K, Satoh N, Satou Y: Gene expression profiles of 
transcription factors and signaling molecules in the ascidian embryo: 
towards a comprehensive understanding of gene networks. Development 
2004, 131:4047-4058. 

60. Choi VM, Harland RM, Khokha MK: Developmental expression of FoxJ1.2, 
FoxJ2, and FoxQI in Xenopus tropicalis. Gene Expr Patterns 2006, 6:443-447. 

61. Santagata S, Resh C, Hejnol A, Martindale MQ, Passamaneck YJ: 
Development of the larval anterior neurogenic domains of Terebratalia 
transversa (Brachiopoda) provides insights into the diversification of 
larval apical organs and the spiralian nervous system. Evodevo 2012, 3:3. 

62. Fritzenwanker JH, Saina M, Technau U: Analysis of forkhead and snail 
expression reveals epithelial-mesenchymal transitions during embryonic 
and larval development of Nematostella vectensis. Dev Biol 2004, 
275:389-402. 

63. Magie CR, Pang K, Martindale MQ: Genomic inventory and expression of 
Sox and Fox genes in the cnidarian Nematostella vectensis. Dev Genes 
Evol 2005,215:618-630. 

64. Momose T, Houliston E: Two oppositely localised frizzled RNAs as axis 
determinants in a cnidarian embryo. PLoS Biol 2007, 5:e70. 

65. Momose T, Derelle R, Houliston E: A maternally localised Wnt ligand 
required for axial patterning in the cnidarian Clytia hemisphaerica. 
Development 2008, 1 35:21 05-21 1 3. 

66. Larroux C, Fahey B, Liubicich D, Hinman VF, Gauthier M, Gongora M, Green 
K, Worheide G, Leys SP, Degnan BM: Developmental expression of 
transcription factor genes in a demosponge: insights into the origin of 
metazoan multicellularity. Evol Dev 2006, 8:150-173. 

67. Castresana J, Feldmaier-Fuchs G, Yokobori S, Satoh N, Paabo S: The mitochondrial 
genome of the hemichordate Balanoglossus carnosus and the evolution of 
deuterostome mitochondria. Genetics 1 998, 1 50:1 1 1 5-11 23. 

68. Bromham LD, Degnan BM: Hemichordates and deuterostome evolution: 
robust molecular phylogenetic support for a hemichordate + echinoderm 
clade. Evol Dev 1 999, 1 :1 66-1 71 . 

69. Cameron CB, Garey JR, Swalla BJ: Evolution of the chordate body plan: 
new insights from phylogenetic analyses of deuterostome phyla. Ptoc 
Natl Acad Sci USA 2000, 97:4469-4474. 

70. Furlong RF, Holland PW: Bayesian phylogenetic analysis supports 
monophyly of ambulacraria and of cyclostomes. Zoolog Sci 2002, 1 9:593-599. 

71. Lowe CJ, Wu M, Salic A, Evans L, Lander E, Stange-Thomann N, Gruber CE, 
Gerhart J, Kirschner M: Anteroposterior patterning in hemichordates and 
the origins of the chordate nervous system. Cell 2003, 1 1 3:853-865. 

72. Aronowicz J, Lowe CJ: Hox gene expression in the hemichordate 
Saccoglossus kowalevskii and the evolution of deuterostome nervous 
systems. Integr Comp Biol 2006, 46:890-901. 

73. Pani AM, Mullarkey EE, Aronowicz J, Assimacopoulos S, Grove EA, Lowe CJ: 
Ancient deuterostome origins of vertebrate brain signalling centres. 
Nature 2012, 483:289-294. 

74. Bateson W: The later stages in the development of Balanoglossus 
kowalevskii with a suggestion as to the affinities of the Enteropneusta. 
0 J Microscrop Sci 1 885, 25:81-1 28. 

75. Schaeffer B: Deuterostome monophyly and phylogeny. Evol Biol 1987, 
21:179-235. 

76. Rychel AL, Smith SE, Shimamoto HT, Swalla BJ: Evolution and development 
of the chordates: collagen and pharyngeal cartilage. Mol Biol Evol 2006, 
23:541-549. 

77. Rychel AL, Swalla BJ: Development and evolution of chordate cartilage. 
J Exp Zool B Mol Dev Evol 2007, 308:325-335. 

78. Okai N, Tagawa K, Humphreys T, Satoh N, Ogasawara M: Characterization 
of gill-specific genes of the acorn worm Ptychodera flava. Dev Dyn 2000, 
217:309-319. 

79. Gillis JA, Fritzenwanker JH, Lowe CG: A stem-deuterostome origin of the 
vertebrate pharyngeal transcriptional network. Proc ft 5oc B 2012, 

279:234-246. 

80. Bullock TH: The anatomical organization of the nervous system of 

Enteropneusta. 0 J Microsc Sci 1 945, 86:55-1 1 1 . 



81 . Knight-Jones EW: On the nervous system of Saccoglossus cambrensis 
(Enteropneusta). Philos Trans R Soc Lond B Biol Sci 1952, 236:315. 

82. Bullock TH: The nervous system of hemichordates. In Structure and 
Function in the Nervous Systems of Invertebrates. San Francisco, CA: WH 
Freeman and Co; 1965:1559-1592. 

83. Lowe CJ, Terasaki M, Wu M, Freeman RM, Runft L, Kwan K, Haigo S, 
Aronowicz J, Lander E, Gruber C, Smith M, Kirschner M, Gerhart J: 
Dorsoventral patterning in hemichordates: insights into early chordate 
evolution. PLoS Biol 2006, 4:e291. 

84. Brown FD, Prendergast A, Swalla BJ: Man is but a worm: chordate origins. 
Genesis 2008, 46:605-613. 

85. Nomaksteinsky M, Rottinger E, Dufour HD, Chettouh Z, Lowe CJ, 
Martindale MQ, Brunet JF: Centralization of the deuterostome nervous 
system predates chordates. Curr Biol 2009, 19:1264-1269. 

86. Kaul S, Stach T: Ontogeny of the collar cord: neurulation in the 
hemichordate Saccoglossus kowalevskii. J Morphol 2010, 271:1240-1259. 

87. Miyamoto N, Wada H: Hemichordate neurulation and the origin of the 
neural tube. Nat Commun 201 3, 4:271 3. 

88. Lowe CJ, Tagawa K, Humphreys T, Kirschner M, Gerhart J: Hemichordate 
embryos: procurement, culture, and basic methods. Methods Cell Biol 
2004, 74:171-194. 

89. Bateson W: The early stages in the development of Balanoglossus (sp. 
incert). 0 J Microscrop Sci 1 884, 24:208-236. 

90. Bateson W: Continued account of the later stages in the development of 
Balanoglossus kowalevskii, with a suggestion as to the affinities of the 
enteropneusta. 0 J Microscrop Sci 1886, 26:511-534. 

91 . Colwin AL, Colwin LH: The normal embryology of Saccoglossus 
kowalevskii. J Morphol 1953, 92:401-453. 

92. Freeman RM, Wu M, Cordonnier-Pratt MM, Pratt LH, Gruber CE, Smith M, 
Lander ES, Stange-Thomann N, Lowe CJ, Gerhart J, Kirschner M: cDNA 
sequences for transcription factors and signaling proteins of the 
hemichordate Saccoglossus kowalevskii: efficacy of the expressed sequence 
tag (EST) approach for evolutionary and developmental studies of a new 
organism. Biol Bull 2008, 214:284-302. 

93. Darras S, Gerhart J, Terasaki M, Kirschner M, Lowe CJ: {betaJ-Catenin 
specifies the endomesoderm and defines the posterior organizer of 
the hemichordate Saccoglossus kowalevskii. Development 201 1, 
138:959-970. 

94. Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, 
Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, Benito-Gutierrez EL, 
Dubchak I, Garcia-Fernandez J, Gibson-Brown JJ, Grigoriev IV, Horton AC, 
de Jong PJ, Jurka J, Kapitonov W, Kohara Y, Kuroki Y, Lindquist E, Lucas S, 
Osoegawa K, Pennacchio LA, Salamov AA, Satou Y, Sauka-Spengler T, 
Schmutz J, Shin-I T, et ah The amphioxus genome and the evolution of 
the chordate karyotype. Nature 2008, 453:1064-1071. 

95. Putnam NH, Srivastava M, Hellsten U, Dirks B, Chapman J, Salamov A, Terry A, 
Shapiro H, Lindquist E, Kapitonov W, Jurka J, Genikhovich G, Grigoriev IV, 
Lucas SM, Steele RE, Finnerty JR, Technau U, Martindale MQ, Rokhsar DS: 
Sea anemone genome reveals ancestral eumetazoan gene repertoire 
and genomic organization. Science 2007, 317:86-94. 

96. Chapman JA, Kirkness EF, Simakov O, Hampson SE, Mitros T, Weinmaier T, 
Rattei T, Balasubramanian PG, Borman J, Busam D, Disbennett K, Pfannkoch C, 
Sumin N, Sutton GG, Viswanathan LD, Walenz B, Goodstein DM, Hellsten U, 
Kawashima T, Prochnik SE, Putnam NH, Shu S, Blumberg B, Dana CE, Gee L, 
Kibler DF, Law L, Lindgens D, Martinez DE, Peng J, et at The dynamic genome 
of Hydra. Nature 201 0, 464:592-596. 

97. Rambaut A: Tree Figure Drawing Tool FigTreei.23. Edinburgh: Institute of 
Evolutionary Biology, University of Edinburgh; 2009. 

98. MrBayes. http://mrbayes.csit.fsu.edu. 

99. Ronquist F, Huelsenbeck JP: MrBayes 3: Bayesian phylogenetic inference 
under mixed models. Bioinformatics 2003, 19:1572-1574. 

1 00. In Le Laboratoire d'lnformatigue, de Robotigue et de Microelectronigue de 
Montpellier (LIRMM). http://phylogeny.lirmm.fr/phylo_cgi/phylogeny.cgi. 

101. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, 
Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O: Phylogeny.fr: robust 
phylogenetic analysis for the non-specialist. Nucleic Acids Res 2008, 
36:W465-W469. 

102. EMBL-EBI. http://www.ebi.ac.uk/Tools/msa/clustalw2/. 

103. Fuchs J, Martindale MQ, Hejnol A: Gene expression in bryozoan larvae 
suggest a fundamental importance of pre-patterned blastemic cells in 
the bryozoan life-cycle. Evodevo 201 1, 2:13. 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



1 04. Katoh M, Katoh M: Human FOX gene family (Review). Int J Oncol 2004, 
25:1495-1500. 

105. Yaklichkin 5, Vekker A, Stayrook S, Lewis M, Kessler DS: Prevalence of the 
EH1 Groucho interaction motif in the metazoan Fox family of 
transcriptional regulators. BMC Genomics 2007, 8:201. 

106. Shimeld SM: A transcriptional modification motif encoded by homeobox 
and fork head genes. FEBS Lett 1997, 410:124-125. 

107. Copley RR: The EH1 motif in metazoan transcription factors. BMC 
Genomics 2005, 6:169. 

108. Colwin AL, Colwin LH: The developmental capacities of separated early 
blastomeres of an enteropneust, Saccoglossus kowalevskii. J Exp Zool 
1950, 115:263-296. 

109. Burdon-Jones C: Development and biology of the larva of 
Saccoglossus horsti (Enteropneusta). Philos Trans R Soc Land B Biol Sci 
1952, 236:553-589. 

110. Rottinger E, Lowe CJ: Evolutionary crossroads in developmental biology: 
hemichordates. Development 2012, 139:2463-2475. 

111. Green SA, Norris RP, Terasaki M, Lowe CJ: FGF signaling induces 
mesoderm in the hemichordate Saccoglossus kowalevskii. Development 
2013, 140:1024-1033. 

112. National Human Genome Research Institute, National Institutes of 
Health, http://research.nhgri.nih.gov/mnemiopsis/. 

113. Marlow H, Tosches MA, Tomer R, Steinmetz PR, Lauri A, Larsson T, Arendt D: 
Larval body patterning and apical organs are conserved in animal 
evolution. BMC Biol 2014, 12:7. 

114. Mango SE, Lambie EJ, Kimble J: The pha-4 gene is required to generate 
the pharyngeal primordium of Caenorhabditis elegans. Development 
1994, 120:3019-3031. 

115. Boyle MJ, Seaver EC: Developmental expression of foxA and gata genes 
during gut formation in the polychaete annelid, Capitella sp I. Evol Dev 
2008, 10:89-105. 

116. Boyle MJ, Seaver EC: Expression of FoxA and GATA transcription factors 
correlates with regionalized gut development in two lophotrochozoan 
marine worms: Chaetopterus (Annelida) and Themiste lageniformis 
(Sipuncula). Evodevo 2010, 1:2. 

1 17. Taguchi S, Tagawa K, Humphreys T, Nishino A, Satoh N, Harada Y: 
Characterization of a hemichordate fork head/HNF-3 gene expression. 
Dev Genes Evol 2000, 210:1 1-17. 

118. Saudemont A, Haillot E, Mekpoh F, Bessodes N, Quirin M, Lapraz F, Duboc V, 
Rottinger E, Range R, Oisel A, Besnardeau L, Wincker P, Lepage T: Ancestral 
regulatory circuits governing ectoderm patterning downstream of Nodal 
and BMP2/4 revealed by gene regulatory network analysis in an 
echinoderm. PLoS Genet 2010, 6:el001259. 

1 19. Friedman JR, Kaestner KH: The Foxa family of transcription factors in 
development and metabolism. Cell Mol Life Sci 2006, 63:231 7-2328. 

1 20. Olsen CL, Jeffery WR: A forkhead gene related to HNF-3beta is required 
for gastrulation and axis formation in the ascidian embryo. Development 
1997, 124:3609-3619. 

121. Martindale MQ, Pang K, Finnerty JR: Investigating the origins of 
triploblasty: 'mesodermal' gene expression in a diploblastic animal, the 
sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). 
Development 2004, 131:2463-2474. 

122. Pohl BS, Knochel W: Of Fox and Frogs: Fox (fork head/winged helix) 
transcription factors in Xenopus development. Gene 2005, 
344:21-32. 

1 23. David ES, Luke NH, Livingston BT: Characterization of a gene encoding a 
developmental^ regulated winged helix transcription factor of the sea 
urchin Strongylocentrotus purpuratus. Gene 1999, 236:97-105. 

124. Minokawa T, Rast JP, Arenas-Mena C, Franco CB, Davidson EH: Expression 
patterns of four different regulatory genes that function during sea 
urchin development. Gene Expr Patterns 2004, 4:449-456. 

1 25. Gamse JT, Sive H: Early anteroposterior division of the presumptive 
neurectoderm in Xenopus. Mech Dev 2001, 104:21-36. 

126. Pohl BS, Knochel S, Dillinger K, Knochel W: Sequence and expression of 
FoxB2 (XFD-5) and Foxllc (XFD-10) in Xenopus embryogenesis. Mech Dev 
2002, 117:283-287. 

127. Pohl BS, Rossner A, Knochel W: The Fox gene family in Xenopus laevis: 
Foxl2, FoxMI and FoxPI in early development. Int J Dev Biol 2005, 
49:53-58. 

128. Mazet F, Shimeld SM: The evolution of chordate neural segmentation. 

Dev Biol 2002, 251:258-270. 



Page 23 of 25 



129. Tan PB, Lackner MR, Kim SK: MAP kinase signaling specificity mediated by 
the LIN-1 Ets/LIN-31 WH transcription factor complex during C. elegans 
vulval induction. Cell 1998, 93:569-580. 

130. Mariani FV, Harland RM: XBF-2 is a transcriptional repressor that converts 
ectoderm into neural tissue. Development 1998, 125:5019-5031. 

131. Seo S, Kume T: Forkhead transcription factors, Foxcl and Foxc2, are 
required for the morphogenesis of the cardiac outflow tract. Dev Biol 
2006, 296:421-436. 

132. Kume T: The cooperative roles of Foxcl and Foxc2 in cardiovascular 
development. Adv Exp Med Biol 2009, 665:63-77. 

1 33. Topczewska JM, Topczewski J, Shostak A, Kume T, Solnica-Krezel L, Hogan BL: 
The winged helix transcription factor Foxcl a is essential for somitogenesis 
in zebrafish. Genes Dev 2001 , 1 5:2483-2493. 

134. Wilm B, James RG, Schultheiss TM, Hogan BL: The forkhead genes, Foxcl 
and Foxc2, regulate paraxial versus intermediate mesoderm cell fate. 
Dev Biol 2004,271:176-189. 

135. Pohl BS, Knochel W: Temporal and spatial expression patterns of FoxD2 
during the early development of Xenopus laevis. Mech Dev 2002, 
111:181-184. 

136. Dirksen ML, Jamrich M: Differential expression of fork head genes during 
early Xenopus and zebrafish development. Dev Genet 1995, 17:107-1 16. 

137. Pohl BS, Knochel W: Overexpression of the transcriptional repressor 
FoxD3 prevents neural crest formation in Xenopus embryos. Mech Dev 

2001, 103:93-106. 

138. Sasai N, Mizuseki K, Sasai Y: Requirement of FoxD3-class signaling for 
neural crest determination in Xenopus. Development 2001, 128:2525-2536. 

139. Lister JA, Cooper C, Nguyen K, Modrell M, Grant K, Raible DW: Zebrafish 
Foxd3 is required for development of a subset of neural crest 
derivatives. Dev Biol 2006, 290:92-104. 

140. Stewart RA, Arduini BL, Berghmans S, George RE, Kanki JP, Henion PD, Look 
AT: Zebrafish foxd3 is selectively required for neural crest specification, 
migration and survival. Dev Biol 2006, 292:174-188. 

141 . Sauka-Spengler T, Bronner-Fraser M: Insights from a sea lamprey into the 
evolution of neural crest gene regulatory network. Biol Bull 2008, 
214:303-314. 

142. Curran K, Lister JA, Kunkel GR, Prendergast A, Parichy DM, Raible DW: 
Interplay between Foxd3 and Mitf regulates cell fate plasticity in the 
zebrafish neural crest. Dev Biol 2010, 344:107-1 18. 

143. Gomez-Skarmeta JL, de la Calle-Mustienes E, Modolell J, Mayor R: Xenopus 
brain factor-2 controls mesoderm, forebrain and neural crest development. 
Mech Dev 1999, 80:15-27. 

144. Odenthal J, Nusslein-Volhard C: Fork head domain genes in zebrafish. Dev 
Genes Evol 1998, 208:245-258. 

145. Lacalli TC, Holland ND, West JE: Landmarks in the anterior central nervous 
system of amphioxus larvae. Philos Trans R Soc Land B Biol Sci 1994, 
344:165. 

146. Holland ND, Panganiban G, Henyey EL, Holland LZ: Sequence and 
developmental expression of AmphiDII, an amphioxus Distal-less gene 
transcribed in the ectoderm, epidermis and nervous system: insights 
into evolution of craniate forebrain and neural crest. Development 1996, 
122:2911-2920. 

147. Lacalli TC: Frontal eye circuitry, rostral sensory pathways and brain 
organization in amphioxus larvae: evidence from 3D reconstructions. 

Philos Trans R Soc Lond B Biol Sci 1 996, 351:243. 

148. Williams NA, Holland PWH: Old head on young shoulders. Nature 1996, 
383:490. 

149. Holland LZ, Holland ND: Evolution of neural crest and placodes: 
amphioxus as a model for the ancestral vertebrate? J Anat 2001, 
199:85-98. 

150. Nash B, Colavita A, Zheng H, Roy PJ, Culotti JG: The Forkhead transcription 
factor UNC-130 is required for the graded spatial expression of the 
UNC-129 TGF-p guidance factor in C. elegans. Genes Dev 2000, 

14:2486-2500. 

151. Sarafi-Reinach TR, Sengupta P: The forkhead domain gene unc-130 
generates chemosensory neuron diversity in C. elegans. Genes Dev 
2000, 14:2472-2485. 

152. Colavita A, Culotti JG: Suppressors of ectopic UNC-5 growth cone steering 
identify eight genes involved in axon guidance in Caenorhabditis elegans. 
Dev Biol 1998, 194:72-85. 

153. Colavita A, Krishna S, Zheng H, Padgett RW, Culotti JG: Pioneer axon 
guidance by UNC-129, a C. elegans TGF-beta. Science 1998, 281:706-709. 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



154. Lapraz F, Besnardeau L, Lepage T: Patterning of the dorsal-ventral axis in 
echinoderms: insights into the evolution of the BMP-chordin signaling 
network. PLoS Biol 2009, 7:el 000248. 

155. Niehrs C: On growth and form: a Cartesian coordinate system of Wnt and 
BMP signaling specifies bilaterian body axes. Development 2010, 137:845-857. 

1 56. Zannini M, Avantaggiato V, Biffali E, Arnone Ml, Sato K, Pischetola M, Taylor 
BA, Phillips SJ, Simeone A, Di Lauro R: TTF-2, a new forkhead protein, 
shows a temporal expression in the developing thyroid which is 
consistent with a role in controlling the onset of differentiation. EMBO J 
1997, 16:3185-3197. 

157. Hiruta J, (vlazet F, Yasui K, Zhang P, Ogasawara M: Comparative expression 
analysis of transcription factor genes in the endostyle of invertebrate 
chordates. Dev Dyn 2005, 233:1031-1037. 

158. Mazet F: The Fox and the thyroid: the amphioxus perspective. Bioessays 

2002, 24:696-699. 

159. El-Hodiri H, Bhatia-Dey N, Kenyon K, Ault K, Dirksen (VI, Jamrich M: Fox 
(forkhead) genes are involved in the dorso-ventral patterning of the 
Xenopus mesoderm. Int J Dev Biol 2001, 45:265-271. 

160. Aitola M, Carlsson P, Mahlapuu M, Enerback S, Pelto-Huikko M: Forkhead 
transcription factor FoxF2 is expressed in mesodermal tissues involved 
in epithelio-mesenchymal interactions. Dev Dyn 2000, 218:136-149. 

161. Mahlapuu M, Ormestad M, Enerback S, Carlsson P: The forkhead 
transcription factor Foxfl is required for differentiation of extra- 
embryonic and lateral plate mesoderm. Development 2001, 128:155-166. 

162. Papalopulu N, Kintner C: A posteriorising factor, retinoic acid, reveals that 
anteroposterior patterning controls the timing of neuronal differentiation 
in Xenopus neuroectoderm. Development 1996, 122:3409-3418. 

163. Hebert JM, Fishell G: The genetics of early telencephalon patterning: 
some assembly required. Nat Rev Neurosci 2008, 9:678-685. 

164. Roth M, Bonev B, Lindsay J, Lea R, Panagiotaki N, Houart C, Papalopulu N: 
FoxGl and TLE2 act cooperatively to regulate ventral telencephalon 
formation. Development 2010, 137:1553-1562. 

165. Ohyama T, Groves AK: Expression of mouse Foxi class genes in early 
craniofacial development. Dev Dyn 2004, 231:640-646. 

166. Stubbs JL, Oishi I, Izpisua Belmonte JC, Kintner C: The forkhead protein 
Foxjl specifies node-like cilia in Xenopus and zebrafish embryos. Nat 
Genet 2008,40:1454-1460. 

167. Blatt EN, Yan XH, Wuerffel MK, Hamilos DL, Brody SL: Forkhead 
transcription factor HFH-4 expression is temporally related to ciliogenesis. 
Am J Respir Cell Mol Biol 1 999, 21:1 68-1 76. 

168. Pohl BS, Knochel W: Isolation and developmental expression of Xenopus 
FoxJl and FoxKl. Dev Genes Evol 2004, 214:200-205. 

169. Yu X, Ng CP, Habacher H, Roy S: Foxjl transcription factors are master 
regulators of the motile ciliogenic program. Nat Genet 2008, 
40:1445-1453. 

1 70. Tichelaar JW, Wert SE, Costa RH, Kimura S, Whitsett JA: HNF-3/forkhead 
homologue-4 (HFH-4) is expressed in ciliated epithelial cells in the 
developing mouse lung. J Histochem Cytochem 1999,47:823-832. 

171. Gomperts BN, Gong-Cooper X, Hackett BP: Foxjl regulates basal body 
anchoring to the cytoskeleton of ciliated pulmonary epithelial cells. 
J Cell So 2004, 1 17:1329-1337. 

172. Nlshimura Y, Hamazaki TS, Komazaki S, Kamimura S, Okochi H, Asashima M: 
Ciliated cells differentiated from mouse embryonic stem cells. Stem Cells 
2006, 24:1381-1388. 

173. Brody SL, Yan XH, Wuerffel (VI K, Song SK, Shapiro SD: Ciliogenesis and 
left-right axis defects in forkhead factor HFH-4-null mice. 71m J Respir 
Cell Mol Biol 2000, 23:45-51. 

174. Zhang M, Bolfing MF, Knowles HJ, Karnes H, Hackett BP: Foxjl regulates 
asymmetric gene expression during left-right axis patterning in mice. 
Biochem Biophys Res Commun 2004, 324:1413-1420. 

1 75. Tamakoshi T, Itakura T, Chandra A, Uezato T, Yang Z, Xue XD, Wang B, 
Hackett BP, Yokoyama T, Miura N: Roles of the Foxjl and Inv genes in the 
left-right determination of internal organs in mice. Biochem Biophys Res 
Commun 2006, 339:932-938. 

176. Tian T, Zhao L, Zhang M, Zhao X, Meng A: Both foxjl a and foxjl b are 
implicated in left-right asymmetric development in zebrafish embryos. 
Biochem Biophys Res Commun 2009, 380:537-542. 

1 77. Chen J, Knowles HJ, Hebert JL, Hackett BP: Mutation of the mouse 
hepatocyte nuclear factor/forkhead homologue 4 gene results in an 
absence of cilia and random left-right asymmetry. J Clin Invest 1998, 
102:1077-1082. 



Page 24 of 25 



178. Sinigaglia C, Busengdal H, Ledere L, Technau U, Rentzsch F: The bilaterian 
head patterning gene six3/6 controls aboral domain development in a 
cnidarian. PLoS Biol 2013, 1 l:el001488. 

1 79. Zhu G, Davis TN: The fork head transcription factor Hcml p participates 
in the regulation of SPC1 1 0, which encodes the calmodulin-binding 
protein in the yeast spindle pole body. Biochim Biophys Acta 1998, 
1448:236-244. 

180. King N, Westbrook MJ, Young SL, Kuo A, Abedln (VI, Chapman J, Fairclough S, 
Hellsten U, Isogai Y, Letunic I: The genome of the choanoflagellate Monosiga 
brevicollis and the origin of metazoans. Nature 2008, 451:783-788. 

181. Nehls M, Pfeifer D, Schorpp M, Hedrich H, Boehm T: New member of the 
winged-helix protein family disrupted in mouse and rat nude mutations. 
Nature 1994, 372:103-107. 

182. Balclunaite G, Keller MP, Balciunaite E, Piali L, Zuklys S, Mathieu YD, Gill J, 
Boyd R, Sussman DJ, Hollander GA: Wnt glycoproteins regulate the 
expression of FoxNI, the gene defective in nude mice. Nat Immunol 
2002,3:1102-1108. 

183. Gouge A, Holt J, Hardy AP, Sowden JC, Smith HK: Foxn4-a new member of 
the forkhead gene family is expressed in the retina. Mech Dev 2001, 
107:203-206. 

184. Li S, Mo Z, Yang X, Price SM, Shen MM, Xiang M: Foxn4 controls the 
genesis of amacrine and horizontal cells by retinal progenitors. Neuron 
2004, 43:795-807. 

185. Kay JN, Baier H: Out-foxing fate; molecular switches create neuronal 
diversity in the retina. Neuron 2004, 43:759-760. 

186. Li S, Misra K, Matise MP, Xiang M: Foxn4 acts synergistically with Mashl to 
specify subtype identity of V2 interneurons in the spinal cord. Proc Natl 
Acad Sci USA 2005, 102:10688-10693. 

1 87. Del Barrio MG, Taveira-Marques R, Muroyama Y, Yuk DI, Li S, Wines-Samuelson M, 
Shen J, Smith HK Xiang M, Rowitch D, Richardson WD: A regulatory network 
involving Foxn4, Mashl and delta-like 4/Notchl generates V2a and 
V2b spinal interneurons from a common progenitor pool. Development 
2007, 134:3427-3436. 

188. Ferland RJ, Cherry TJ, Preware PO, Morrisey EE, Walsh CA: Characterization 
of Foxp2 and Foxpl mRNA and protein in the developing and mature 
brain. J Comp Neurol 2003, 460:266-279. 

189. Bonkowsky JL, Chien CB: Molecular cloning and developmental 
expression of foxP2 in zebrafish. Dev Dyn 2005, 234:740-746. 

190. Schon C, Wochnik A, Rossner A, Donow C, Knochel W: The FoxP subclass 
in Xenopus laevis development. Dev Genes Evol 2006, 216:641-646. 

191 . Tamura S, Morikawa Y, Miyajima A, Senba E: Expression of oncostatin M in 
hematopoietic organs. Dev Dyn 2002, 225:327-331. 

192. Shen X, Cui J, Nagahama Y: The forkhead gene family in medaka: 
expression patterns and gene evolution. Cytogenet Genome Res 2012, 
136:123-130. 

193. Lai CS, Fisher SE, Hurst JA, Vargha-Khadem F, Monaco AP: A forkhead- domain 
gene is mutated in a severe speech and language disorder. Nature 2001, 
413:519-523. 

194. Fisher SE, Scharff C: FOXP2 as a molecular window into speech and 
language. Trends Genet 2009, 25:166-177. 

195. Takahashi H, Takahashi K, Liu FC: FOXP genes, neural development, 
speech and language disorders. Adv Exp Med Biol 2009, 665:1 17-129. 

196. Hannenhalli S, Kaestner KH: The evolution of Fox genes and their role in 
development and disease. Nat Rev Genet 2009, 10:233-240. 

197. Santos ME, Athanasiadis A, Leitao AB, DuPasquier L, Sucena E: Alternative 
splicing and gene duplication in the evolution of the FoxP gene 
subfamily. Mol Biol Evol 201 1, 28:237-247. 

198. Rottinger E, Martindale MQ: Ventralization of an indirect developing 
hemichordate by NiCI 2 suggests a conserved mechanism of dorso-ventral 
(D/V) patterning in Ambulacraria (hemichordates and echinoderms). Dev 
Biol 2011,354:173-190. 

199. Nielsen C: Animal evolution: interrelationships of the living phyla. New York: 
Oxford University Press; 1995. 

200. Jagersten G: The Evolution of the Metazoan Life Cycle, A comprehensive 
Theory. Waltham, MA: Academic Press; 1972. 

201. Peterson KJ, Cameron RA, Davidson EH: Set-aside cells in maximal indirect 
development: evolutionary and developmental significance. Bioessays 
1997, 19:623-631. 

202. Rentzsch F, Fritzenwanker JH, Scholz CB, Technau U: FGF signalling 
controls formation of the apical sensory organ in the cnidarian 
Nematostella vectensis. Development 2008, 135:1761-1769. 



Fritzenwanker et at. EvoDevo 2014, 5:17 
http://www.evodevojournal.eom/content/5/1/17 



Page 25 of 25 



203. Nielsen C: Animal evolution: interrelationships of the living phyla. 3rd edition. 
New York: Oxford University Press; 2012. 

204. Wolpert L: From egg to adult to larva. Evol Dev 1 999, 1 :3-4. 

205. Raff RA: Origins of the other metazoan body plans: the evolution of 
larval forms. Philos Trans R Soc Land B Biol Sci 2008, 363:1473-1479. 

206. Yaguchi S, Yaguchi J, Wei Z, Shiba K, Angerer LM, Inaba K: ankAT-1 is a 
novel gene mediating the apical tuft formation in the sea urchin 
embryo. Dev Biol 2010, 348:67-75. 

207. Satir P, Christensen ST: Overview of structure and function of mammalian 
cilia. Annu Rev Physiol 2007, 69:377-400. 

208. Sasai N, Briscoe J: Primary cilia and graded Sonic Hedgehog signaling. 
Wiley Interdiscip Rev Dev Biol 201 2, 1 :753-772. 

209. Thomas J, Morle L, Soulavie F, Laurencon A, Sagnol S, Durand B: 
Transcriptional control of genes involved in ciliogenesis: a first step in 
making cilia. Biol Cell 2010, 102:499-513. 

210. Vij S, Rink JC, Ho HK, Babu D, Eitel (VI, Narasimhan V, Tiku V, Westbrook J, 
Schierwater B, Roy S: Evolutionarily ancient association of the FoxJI 
transcription factor with the motile ciliogenic program. PLoS Genet 2012, 
8:el003019. 

21 1. Wotton KR, Shimeld SM: Analysis of lamprey clustered Fox genes: insight 
into Fox gene evolution and expression in vertebrates. Gene 2011, 

489:30-40. 

212. Dominguez P, Jacobson AG, Jefferies RP: Paired gill slits in a fossil with a 
calcite skeleton. Nature 2002, 417:841-844. 

213. Ogasawara M, Wada H, Peters H, Satoh N: Developmental expression of 
Paxl/9 genes in urochordate and hemichordate gills: insight into 
function and evolution of the pharyngeal epithelium. Development 1999, 
126:2539-2550. 

214. Nissen RM, Yan J, Amsterdam A, Hopkins N, Burgess SM: Zebrafish foxi one 
modulates cellular responses to Fgf signaling required for the integrity 
of ear and jaw patterning. Development 2003, 1 30:2543-2554. 

215. Solomon KS, Kudoh T, Dawid IB, Fritz A: Zebrafish foxil mediates otic 
placode formation and jaw development. Development 2003, 1 30:929-940. 

216. De Robertis EM, Sasai Y: A common plan for dorsoventral patterning in 
Bilateria. Nature 1 996, 380:37-40. 

217. Arendt D, Nubler-Jung K: Dorsal or ventral: similarities in fate maps and 
gastrulation patterns in annelids, arthropods and chordates. Mech Dev 
1997, 61:7-21. 

218. Gerhart J: Inversion of the chordate body axis: are there alternatives? 

Proc Natl Acad Sci USA 2000, 97:4445-4448. 

219. Gerhart J: The deuterostome ancestor. J Cell Physiol 2006, 209:677-685. 

220. Steiner AB, Engleka MJ, Lu Q, Piwarzyk EC, Yaklichkin S, Lefebvre JL, 
Walters JW, Pineda-Salgado L, Labosky PA, Kessler DS: FoxD3 regulation 
of Nodal in the Spemann organizer is essential for Xenopus dorsal 
mesoderm development. Development 2006, 133:4827-4838. 

221. Kozmikova I, Smolikova J, Vlcek C, Kozmik Z: Conservation and 
diversification of an ancestral chordate gene regulatory network for 
dorsoventral patterning. PLoS One 2011, 6:el4650. 



doi:10.1 186/2041-9139-5-17 

Cite this article as: Fritzenwanker et ai: The Fox/Forkhead transcription 
factor family of the hemichordate Saccoglossus kowalevskii. EvoDevo 

2014 5:17. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f \ n! _, ul _-| r Q r,tr=l 

www.biomedcentral.com/submit B,onne «« central 



V 



